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1  Dispersion  relation  for  a  cold  two  ion  component  plasma  for 
nonparallel  propagation.  The  various  branches  are  labeled  ac¬ 
cording  to  their  predominate  polarization  with  RHCP  dashed 
and  LHCP  solid.  Coupling  occurs  at  the  crossover  frequency 
(u>cr)  where  the  polarization  of  branches  I  and  II  switches 
abruptly.  Tunneling  occurs  across  the  shaded  “gap”  between 
the  cutoff  frequency  (uco)  and  the  heavy  ion  gyrofrequency 

(ftd) . 

2  The  locus  of  solutions  to  Eq.  (1.11)  for  varying  k.  The  co¬ 
ordinate  z  corresponds  to  the  ratio  For  k  =  0,  the 

condition  is  satisfied  at  crossover  frequency  and  at  zero  fre¬ 
quency  while  for  large  values  of  k  the  condition  is  satisfied 
near  the  oxygen  resonances  and  the  hydrogen  resonance  (not 
shown).  The  arrows  indicate  the  direction  of  increasing  k. 
For  very  small  densities  (less  than  0.2%)  the  topology  changes 
and  the  crossover  frequency  is  connected  to  the  oxygen  gy¬ 
rofrequency  whereas  the  zero  frequency  solution  is  connected 
to  the  hydrogen  gyrofrequency . 

3  Illustration  of  the  meaning  of  the  coefficients  for  a  downgoing 

ion  cyclotron  wave.  Waves  are  reflected  back  along  the  two 
upward  propagating  modes  while  a  wave  is  transmitted  on  the 
downgoing  magnetosonic  branch . 

4  Illustration  of  the  meaning  of  the  coefficients  for  a  downgoing 

magnetosonic  wave.  Waves  are  reflected  back  along  the  two 
upward  propagating  modes  while  a  wave  is  transmitted  on  the 
downgoing  magnetosonic  branch . 

5  Sectors  of  the  complex  plane  divided  according  to  regions  of 

dominance  and  subdominance  of  the  asymptotic  WKB  solu¬ 
tions.  The  Stokes  lines  and  anti-Stokes  lines  are  solid  and 
dashed  respectively . 

6  Illustration  of  the  physical  meaning  of  the  coefficients  for  a 

downgoing  ion-cyclotron  wave  for  large  angle  of  propagation. 
Waves  are  transmitted  onto  mode  III  and  reflected  back  along 
mode  II . 
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the  real  axis  to  some  point  p.  To  continue  the  solution  to  the 
positive  real  axis,  we  take  a  contour  along  a  circular  arc  above 
the  pole.  Finally  the  solution  may  be  obtained  by  continuing 


the  solution  along  the  real  axis . 40 

8  In  regions  far  from  the  pole  the  WKB  solutions  are  asymp¬ 
totically  valid.  The  numerical  solution  connects  the  regions  in 
which  the  WKB  solutions  are  accurate . 41 


9  For  the  limit  of  large  “gap”  the  numerical  solutions  and  the 
analytical  results  are  in  reasonable  agreement.  We  have  taken 
e=0.03  and  ij  =  5%  for  which  there  is  no  absorption.  The 
numerical  solutions  are  plotted  in  heavy  set  type  and  the  an¬ 
alytical  results  are  indicated  in  light  type.  Note  that  for  the 
numerical  results  the  coefficient  exp(— nfi+fc)  although  small 
is  nonzero  in  contrast  to  the  analytical  results.  This  behavior 
is  the  result  of  coupling  between  the  two  Alfven  waves  below 
the  ion  resonance  corresponding  to  the  coupling  point  at  zero 
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Part  1 

An  Analytical  Study  of  Mode 
Conversion 

1  Introduction 

It  is  well  known  that  wave  propagation  in  a  multiple  species  plasma  is  sub¬ 
stantially  modified  near  the  ion  reson&nces[Buchsbaum ,  1960a,  1960b].  This 
modification  leads  to  important  effects  which  can  be  utilized,  for  example, 
as  a  diagnostic  to  determine  plasma  concentration[Gwrne«  et  al.,  1965]  or  as 
a  tool  to  heat  plasma  via  ion-cyclotron  resonance  heating  (ICRH)[5wanson, 
1976;  Jacquinot  et  al.,  1977;  Perkins,  1977].  The  effects  of  a  minority  plasma 
are  substantial  in  that  modes  which  would  not  ordinarily  heat  cold  ions  may 
energize  the  minority  population.  In  this  report,  we  address  wave  propaga¬ 
tion  in  a  multispecies  plasma  in  the  context  of  a  specific  problem  of  impor¬ 
tance  in  space  physics.  In  this  problem,  waves  propagate  from  a  region  of 
small  magnetic  field  into  a  region  of  increasing  magnetic  field.  In  some  local¬ 
ized  region  the  waves  encounter  the  local  minority  ion  gyrofrequency  (oxygen 
in  this  case)  and  wave  power  is  reflected,  transmitted,  mode  converted,  and 
absorbed.  Under  reasonable  conditions,  the  wave  power  absorbed  by  the 
minority  species  is  sufficient  to  energize  ions  to  observed  energies. 

One  of  the  more  intriguing  phenomena  of  space  physics  is  the  energization 
of  ionospheric  plasma  to  magnetospheric  energies.  A  particular  manifestation 
of  this  outflow  are  the  events  termed  ion  conics.  Ion  conics  are  essentially 
the  signature  of  a  heating  process  transverse  to  the  magnetic  field.  The 
conic  shape  (distribution  with  a  fixed  pitch  angle  in  velocity  space)  may 
be  heuristically  understood  in  terms  of  conservation  of  the  first  adiabatic 
invariant  where  perpendicular  energy  is  converted  into  parallel  energy  in  the 
geometry  of  a  decreasing  magnetic  field. 

One  particular  variety  of  conic  which  has  provided  remarkable  confir¬ 
mation  of  theory  is  that  which  is  observed  along  closed  field  lines  which 
emanate  from  the  auroral  region  into  the  central  plasma  sheet  (CPS).  These 
conics  are  generally  observed  in  association  with  a  broadband  spectrum  of 
downcoming  waves  which  peaks  at  low  frequencies  and  contains  substantial 


wave  energies  at  the  local  oxygen  cyclotron  frequency  [ Gurnett  et  al. ,  1984]. 
Chang  et  al.  [1986]  suggested  a  heuristic  model  in  which  these  waves  interact 
resonantly  with  ions  over  an  extended  altitude  corresponding  to  a  range  of 
cyclotron  resonance  with  the  broad  range  of  frequencies  within  the  observed 
spectrum.  The  energies  and  pitch  angles  of  the  conics  were  in  agreement 
with  this  model.  Moreover,  a  more  rigorous  analysis  [Retterer  et  <z/.,  1987; 
Crew  and  Chang ,  1988;  Crew  et  a/.,  1990]  not  only  confirmed  this  theory,  but 
provided  what  may  be  the  best  agreement  between  theory  and  observation  in 
space  physics.  Perhaps  the  most  remarkable  aspect  of  this  analysis  is  that  it 
contained  only  one  free  parameter,  namely  the  fraction  of  left-hand  circular 
polarization  (LHCP)  (i.e.  the  amount  of  wave  power  which  can  resonate  with 
the  ions).  Although  the  mechanism  by  which  the  conics  are  formed  seems  to 
be  well  understood,  the  origin  and  nature  of  the  ambient  turbulence  remain 
somewhat  uncertain. 

The  existence  of  these  waves  has  remained  somewhat  of  a  mystery  in 
that  careful  scrutiny  of  the  data  has  to  date  precluded  a  local  source  (very 
little  free  energy  is  available  in  this  region  of  the  magnetosphere).  In  light 
of  the  observations  associated  with  the  conic  events,  we  suggested  a  scenario 
[Johnson  et  al. ,  1989]  in  which  the  waves  are  generated  in  the  equatorial  re¬ 
gion  of  the  magnetosphere,  propagate  along  magnetic  fields  to  lower  altitudes 
(larger  magnetic  field)  and  are  then  absorbed,  to  some  extent,  at  the  local 
cyclotron  resonance.  In  this  manner,  the  minority  oxygen  is  energized  to 
magnetospheric  energies.  According  to  this  scenario  as  well  as  direct  satel¬ 
lite  measurement  of  the  Poynting  flux,  the  source  of  the  observed  spectrum 
would  appear  to  be  from  “above”  the  region  in  which  ion-heating  occurs. 
That  is,  the  waves  propagate  from  regions  of  small  magnetic  field  into  re¬ 
gions  of  large  magnetic  field.  An  obvious  question  then  arises.  How  much 
of  the  wave  energy  is  absorbed  by  ions  at  the  local  oxygen  resonance,  and  is 
that  absorption  sufficient  to  energize  the  ions  to  magnetospheric  energies? 

The  purpose  of  this  present  work  is,  therefore,  to  demonstrate  that  the 
assumed  absorption  of  wave  energy  by  the  ion  conics  is  appreciable  within 
the  framework  of  a  reasonable  model.  In  order  to  determine  the  amount  of 
wave  absorption,  we  simplify  the  geometry  and  make  suitable  approximations 
which  retain  the  most  essential  gradients  involved  in  the  wave  equations. 

It  will  be  obvious  from  the  structure  of  these  equations  that  the  four  low 
frequency  cold  plasma  modes  couple  strongly  near  the  minority  ion  gyrofre- 
quency.  Hence,  power  may  be  exchanged  between  the  modes  so  that  both  the 
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Figure  1:  Dispersion  relation  for  a  cold  two  ion  component  plasma  for 
nonparallel  propagation.  The  various  branches  are  labeled  according  to  their 
predominate  polarization  with  RHCP  dashed  and  LHCP  solid.  Coupling 
occurs  at  the  crossover  frequency  (t^Cr)  where  the  polarization  of  branches  I 
and  II  switches  abruptly.  Tunneling  occurs  across  the  shaded  “gap”  between 
the  cutoff  frequency  ( loco )  and  the  heavy  ion  gyrofrequency  (flci). 

magnetosonic  and  shear  Alfven  modes  can  undergo  substantial  absorption. 
We  analyze  these  mode  conversion  problems  using  a  variety  of  techniques. 
We  obtain  analytical  results  using  a  perturbation  analysis,  the  saddle-point 
method,  and  the  phase-integral  method,  and  we  discuss  parallel  work  which 
employs  an  asymptotic  expansion  of  the  Fourier  transform  of  the  differential 
equation.  The  various  approaches  each  provide  some  insight  as  to  the  physics 
described  by  the  equations.  In  part  2  of  this  report  we  present  a  numerical 
solution  which  provides  not  only  a  check  of  the  analytical  results,  but  also 
serves  to  retain  the  presence  of  multiple  saddle  points  which  are,  in  some 
cases,  not  well  approximated  by  an  analytical  solution. 

In  order  to  understand  more  completely  the  processes  which  we  are  exam¬ 
ining,  consider  the  dispersion  relation  of  Figure  1.  This  dispersion  relation 
is  characterized  by  three  separate  modes  which  we  have  labeled  I,  II,  and 
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III.  Mode  I  is  the  magnetosonic  branch  and  is  predominately  right-hand  cir¬ 
cularly  polarized  at  high  frequency  whereas  mode  II,  which  is  the  hydrogen 
ion-cyclotron  branch,  is  predominately  LHCP.  Below  the  oxygen  cyclotron 
frequency  mode  II  continues  to  propagate  as  the  right-hand  circularly  po¬ 
larized  (RHCP)  magnetosonic  branch  and  mode  III  represents  the  LHCP 
ion-cyclotron  branch.  Near  the  crossover  frequency  (the  relevant  frequen¬ 
cies  in  a  multispecies  plasma  are  discussed  in  Appendix  C),  modes  I  and  II 
couple  and  exchange  energy  [Smith  and  Brice ,  1964],  Mode  I  also  suffers  a 
cutoff  at  the  cutoff  frequency  and  is  separated  from  mode  III  by  an  evanes¬ 
cent  region.  As  waves  move  into  regions  of  larger  magnetic  field  they  move 
“down”  along  the  dispersion  curves  according  to  the  WKB  approximation  in 
that  they  satisfy  the  local  dispersion  relation.  The  scenario  in  which  waves 
propagate  earthward  into  a  region  of  larger  magnetic  field  (as  is  the  case  for 
waves  propagating  from  the  equatorial  region  to  the  auroral  region)  corre¬ 
sponds  to  waves  incident  from  “above”  the  coupling-evanescent  region.  Some 
of  the  incident  wave  power  is  then  reflected  back  along  branches  I  and  II, 
some  wave  power  is  transmitted  along  branches  II  and  III  “below”  the  heavy 
ion  cyclotron  frequency,  and  some  wave  power  is  absorbed  in  the  evanescent 
region.  Particle  heating  may  result  from  ion-cyclotron  damping  of  branch  III 
(for  a  thermal  plasma)  and  from  the  calculated  absorption  which  typically 
corresponds  to  mode  conversion  to  a  thermal  plasma  wave  which  is  strongly 
damped. 


2  Basic  Equations 

Our  goal  is  to  describe  the  behavior  of  low  frequency  waves  near  the  series 
of  crossover  cutoff  and  resonance  frequencies  found  in  a  multispecies  plasma. 
In  the  CPS  the  typical  Alfven  wavelength  far  exceeds  even  the  most  ener¬ 
getic  ion  gyroradius  so  that  a  cold  plasma  description  is  reasonable.  The 
essential  feature  that  we  wish  to  retain  is  the  variation  of  the  magnetic  field 
strength  along  field  lines  (density  gradients  do  not  affect  the  location  of  the 
resonances  and  only  density  ratios  affect  the  location  of  the  other  special 
frequencies).  The  differential  equations  are  most  tractable  and  still  retain 
the  essential  features  which  we  wish  to  incorporate  into  our  analysis  if  we 
consider  a  plasma  in  which  the  mode  conversion  processes  are  dominated  by 
the  variation  in  the  strength  of  the  magnetic  field  along  one  dimension  and 


4 


take  the  density  to  be  constant. 

At  low  altitudes  where  the  heating  and  mode  conversion  processes  are 
presumed  to  occur,  the  primary  variation  in  the  magnetic  field  strength  is 
along  the  direction  of  the  magnetic  field.  Thus,  we  assume  that  the  only 
magnetic  field  variations  are  along  the  field  line,  and  we  neglect  the  variations 
perpendicular  to  the  field  line  and,  in  effect,  the  field  line  curvature.  Such 
an  approximation  is  reasonable  if  the  conversion  process  is  localized  in  some 
sense  with  respect  to  the  scale  length  of  the  magnetic  field.  Indeed,  the  scale 
over  which  the  mode  conversion  occurs  is  asymptotically  small  compared 
with  variations  in  the  magnetic  field  so  that  such  an  approximation  is  valid. 

We  should  emphasize  that  perpendicular  gradients  are  weak  compared 
to  parallel  gradients  for  this  problem  so  that  their  effects  are  not  important 
in  our  calculations.  Thus,  our  work  is  different  from  other  studies  in  which 
perpendicular  gradients  dominate  [ Swanson ,  1976;  Jacquinot  et  al .,  1977; 
Perkins,  1977].  We  should  point  out  that  variations  perpendicular  to  the 
magnetic  field  do  not  couple  the  two  propagating  cold  plasma  modes  in  that 
wave  propagation  across  the  field  is  described  by  a  second  order  differential 
equation.  However,  it  should  be  noted  that  variations  across  the  field  should 
give  rise  to  coupling  between  cold  plasma  waves  and  thermal  ion-ion  harmonic 
waves  [Sicanson,  1976;  Jacquinot  et  al.,  1977;  Perkins,  1977]. 

Maxwell’s  equations  combined  with  the  momentum  transfer  equations 
describe  wave  propagation  in  a  plasma  medium  [ Stix ,  1962].  In  light  of 
the  approximations  detailed  above,  these  equations  reduce  to  a  set  of  cou¬ 
pled  ordinary  differential  equations  for  the  components  of  the  electric  field 
[Forsterling,  1942].  For  low  frequencies,  the  electric  field  parallel  to  the  mag¬ 
netic  field  becomes  negligible  [Swanson,  1989]  so  that  the  equations  for  the 
circularly  polarized  electric  field  components  perpendicular  to  the  magnetic 
field,  E±  =  Ex  ±  iEy  take  the  compact  form 

e2tf"  +  M'k  =  0  (1.1) 

where 

M 

The  vector  $  consists  of  two  components  E _  and  E+  which  we  will  also 
refer  to  as  0  and  0.  The  characteristics  of  the  medium  are  contained  in  the 
parameters:  e  =  1  /k^Ls  ,  «  =  k^/kA,  kA  —  oj/va,  Lb  —  B~ldB / dz.  The 


r-K2/ 2  k2/2 

k2/ 2  1-k2/  2 


(1.2) 
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rescaled  Stix  functions 


r  _  (a?±a>c0) 

/  cuC0  (u>±Qci)(u>±ftc2)  (  ' 

where  is  the  Alfven  velocity.  We  have  taken  the  Fourier  transform  in  the 
direction  perpendicular  to  the  magnetic  field.  In  addition,  we  have  rescaled 
the  coordinate  variable  to  some  scale  length,  Lg ,  which  appropriately  de¬ 
scribes  the  scale  length  of  the  magnetic  field. 

These  equations  are  related  to  a  dispersion  relation  which  is  obtained 
from  the  Fourier  transform  of  Eq.  (1.1).  Indeed,  the  substitution  itd/dz  n 
recovers  the  related  dispersion  relations 


n± 


2 


r  +  / 
2 


K 


~Y±  ( 


^  /  2 

— ) 

Q  / 


(1.4) 


We  can  make  a  number  of  observations  concerning  the  differential  Eq. 
(1.1).  These  equations  are  primarily  characterized  by  the  quantities  e,  k, 
and  the  relative  spacing  of  the  cyclotron  frequencies  and  the  cutoff.  The 
essence  of  the  WI<B  approximation  is  contained  in  the  parameter  e.  For 
small  e  the  WKB  solutions  most  nearly  approximate  the  correct  solution  in 
the  regions  where  r  and  /  are  well  behaved.  The  parameter  k,  on  the  other 
hand  characterizes  the  coupling  between  the  two  modes.  Indeed,  from  Eq. 
(1.1,  1.2)  it  is  obvious  that  for  k  —  0,  the  two  equations  are  completely 
uncoupled  and  each  mode  may  be  described  separately.  In  this  case,  the 
equations  simply  describe  the  parallel  propagating  R  and  L  modes. 


3  Analytical  Considerations 

Mode  conversion  is  generally  characterized  by  the  coalescence  of  two  propa¬ 
gating  modes  [ Budden ,  1985;  Fuchs  et  al. ,  1981;  Swanson ,  1989].  In  a  region 
where  mode  conversion  occurs,  the  propagating  modes  which  coalesce  may, 
in  fact,  exchange  substantial  energy.  The  local  properties  of  a  medium  in  one 
dimension  are  characterized  by  the  dispersion  relation  D(k,z)  =  0.  Saddle 
points  of  the  dispersion  relation  where  at  least  two  roots  coalesce  occur  where 
dD/dk  =  0  which  is  satisfied  for  the  mapping  k  ~  kc(z)  which  depends  upon 
the  spatial  coordinate,  z.  At  a  point,  z&,  where  the  dispersion  relation  is 
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satisfied  kc(zb )  becomes  a  double  root  of  the  dispersion  relation.  The  point, 
Zb,  is  referred  to  as  a  branch  point.  Near  the  coupling  point,  the  dispersion 
equation  may  be  described  in  terms  of  an  embedded  dispersion  relation 

D{kc,  z)  +  -  kc)2Dkk{kc,z)  =  0  (1.5) 

Then  a  differential  equation  describing  the  two  mode  coupling  may  be  ob¬ 
tained  from  the  Fourier  transform  in  k  so  that 


rt>'\z)  +  Q(z)if>  =  0 


(1.6) 


with  Q(z )  =  —2 D(kc,  z)/ Dkk(kc,  z).  By  definition,  the  potential  Q  vanishes 
at  branch  points  Zb  where  D  =  0.  When  coupling  points  are  well  separated, 
the  embedded  equation  contains  sufficient  information  to  describe  the  cou¬ 
pling  between  the  embedded  modes.  Such  coupling  is  typically  very  strong  if 
the  coupling  point  Zb  occurs  close  to  the  real  axis  and  is  weak  if  the  coupling 
point  is  far  from  the  real  axis. 

The  particular  dispersion  relation  that  we  are  considering  for  low  fre¬ 
quency  waves  may  be  obtained  from  the  Fourier  transform  of  Eq.  (1.1) 

D(k,  z)  =  ( k 2  -  r)(k2  -/)  +  K2(k2  -  r-±±)  (1.7) 

where  it  is  understood  that  r  and  /  are  explicit  functions  of  z.  The  saddle 
points  are  located  along  contours 

kc{z)  =  0,  ^(r  +  l  -  k2)  (1.8) 

The  branch  points  in  space  related  to  these  saddle  points  may  be  obtained  by 
solving  the  dispersion  relation  (1.7)  along  the  saddle  point  contours,  kc(z). 
For  the  coupling  point  at  kc  =  0  we  find  the  cutoff  condition 


r'  =  y(r  +  0  (1.9) 

which  occurs  at  the  position  at  which  the  magnetic  field  satisfies  the  condition 

—  k2/ 2 

it ,  ,2  „2,  .2  /O  (I.IUJ 


U) 


2  2 
u?  or 


ul  -  2 
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This  is  simply  the  cutoff  condition  for  the  dispersion  relation.  For  small 
values  of  k,  this  cutoff  occurs  at  the  cutoff  frequency,  u>co,  which  for  parallel 
propagation  is  the  cutoff  of  the  LHCP  L-mode.  Near  k2  =  2  the  denominator 
diverges.  (Actually  the  root  approaches  the  lower  hybrid  frequency,  but  our 
approximations,  suitable  for  low  frequencies  do  not  adequately  describe  this.) 
On  the  other  hand,  for  k  >  y/2  a  cutoff  appears  at  low  frequencies  which 
increases  as  a  function  of  k  up  to  the  ion-ion  resonance.  At  some  critical  k 
the  cutoff  passes  through  the  resonance. 

The  coupling  point  kc  =  (r  +  /  —  k2)/2,  when  substituted  into  the  disper¬ 
sion  relation,  yields  the  condition 

d2  +  ^=0  (l.n) 

where  d  is  the  difference  function 

J  =  r  ~~  [  _  ffclf^c2  Lo{iQ2  —  Lp'lr) 

2  ujco  1  J 

This  condition  simply  corresponds  to  a  double  root  of  the  dispersion  relation 
for  the  modes  given  in  Eq.  (1.4).  For  k  —  0,  this  condition  is  satisfied  at  the 
cutoff  and  crossover  frequencies  which  are  the  zeros  of  the  function  d.  As 
K  oo,  this  condition  may  only  be  satisfied  near  the  ion  gyrofrequencies. 
For  ft  7^  0  it  is  clear  that  this  condition  may  be  satisfied  only  for  values 
of  complex  d  which  correspond  to  complex  values  of  the  magnetic  field  in 
which  case  the  branch  points,  z\>,  are  complex.  It  should  also  be  clear  from 
the  complex  conjugate  of  Eq.  (1.11)  that  the  complex  branch  points  always 
occur  in  complex  conjugate  pairs.  In  Figure  2  we  plot  values  of  the  complex 
frequency  along  which  Eq.  (1.11)  is  satisfied.  Eq.  (1.11)  is  an  eighth  order 
equation  in  frequency  with  four  complex  conjugate  pair  solutions  except  at 
k  =  0  and  ft  ^  oo  where  the  modes  coalesce  into  a  pair  of  real  double 
roots.  For  reasonable  ionospheric  values  of  the  oxygen  density,  the  crossover 
frequency  connects  to  the  oxygen  cyclotron  frequency  and  the  double  root 
at  u  =  0  connects  to  the  heavy  ion  gyrofrequency.  From  the  topology,  we 
expect  coupling  to  be  largest  when  the  coupling  points  lie  near  the  real  axis. 
In  this  sense,  we  see  that  the  largest  coupling  takes  place  for  small  values  of 
ft  near  the  crossover  frequency  and  near  the  u  =  0  root  (which  corresponds 
physically  to  very  low  altitudes).  On  the  other  hand,  for  large  values  of  ft 
the  coupling  occurs  near  the  two  gyrofrequencies  (the  other  gyrofrequency 
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Re(z) 


Figure  2:  The  locus  of  solutions  to  Eq.  (1.11)  for  varying  k.  The  coordinate 
z  corresponds  to  the  ratio  a>/flci.  For  k  =  0,  the  condition  is  satisfied  at 
crossover  frequency  and  at  zero  frequency  while  for  large  values  of  k  the 
condition  is  satisfied  near  the  oxygen  resonances  and  the  hydrogen  resonance 
(not  shown).  The  arrows  indicate  the  direction  of  increasing  «.  For  very  small 
densities  (less  than  0.2%)  the  topology  changes  and  the  crossover  frequency 
is  connected  to  the  oxygen  gyrofrequency  whereas  the  zero  frequency  solution 
is  connected  to  the  hydrogen  gyrofrequency. 


at  z  =  16  is  not  shown  in  Figure  1.2).  The  topology  of  these  branch  points 
is  discussed  in  detail  by  Le  Queau  and  Roux  [1992], 

In  the  following  sections  we  will  carefully  consider  the  problem  at  hand, 
and  using  suitable  approximations,  we  will  extract  the  behavior  of  the  cou¬ 
pling  coefficients  analytically.  In  part  2  of  this  report,  we  solve  the  differential 
Eq.  (1.1)  numerically  and  obtain  coefficients  valid  for  all  relevant  k. 

3.1  Mode  Coupling  at  Crossover  Frequency 

For  k<1,  the  two  modes  coalesce  near  the  crossover  frequency  and  couple 
strongly.  We  may  describe  this  coupling  in  a  number  of  ways.  In  doing  so, 
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we  shall  demonstrate  that  a  substantial  amount  of  energy  will  be  exchanged 
between  the  two  propagating  modes  near  the  crossover  frequency  for  k  <C  1. 
For  k  >  1  the  magnetosonic  mode  is  cut  off  at  the  lower  hybrid  frequency 
which  lies  above  the  light  ion  cyclotron  frequency  so  that  only  two  modes 
propagate  near  the  heavy  ion  gyrofrequency.  The  character  of  the  problem 
therefore  changes  as  the  value  of  k  changes  across  some  critical  value.  In  this 
section,  we  will  concentrate  on  values  of  k  below  this  critical  value  such  that 
four  modes  propagate  in  the  vicinity  of  the  heavy  ion  gyroresonance.  In  a 
later  section,  we  address  propagation  for  larger  values  of  k. 

One  way  in  which  to  understand  coupling  near  the  crossover  frequency  is 
to  consider  a  perturbation  about  parallel  propagation  [Johnson  et  al.,  1989; 
Budden,  1985].  Such  an  analysis  is  particularly  insightful  because  it  is  then 
quite  clear  that  the  coupling  results  from  the  presence  of  finite  k.  Indeed, 
the  coupling  between  the  two  modes  is  found  from  Eq.  (1.1)  to  be  of  order 

AC2. 

For  small  k  the  fields  may  be  expanded 

E±  =  E^  +  k2E^  +  ...  (1.13) 

Eq.  (1.1)  then  reduces  to  a  hierarchy  of  second  order  differential  equations 
which  are  coupled  by  the  known  lower  order  fields.  The  zeroth  order  so¬ 
lutions  are  the  two  uncoupled  parallel  propagating  R,L-modes.  The  higher 
order  fields  may  be  determined  using  variation  of  parameters.  In  order  to 
determine  the  behavior  of  a  downcoming  ion-cyclotron  (magnetosonic)  wave, 
we  must  prescribe  that  E 1°^  (E+^)  vanishes.  In  either  case  the  first  order 
solution  for  the  coupled  magnetosonic  (ion-cyclotron)  wave  is  found  to  have 
the  amplitude 

A(z)  =  i h  />['(^)r^p  (-)/(< w?)]1/2  -  [((?)],/2)V)  (i.i4) 

where  h  is  some  suitably  chosen  reference  altitude  above  the  coupling  region. 
The  amplitude  of  the  coupled  wave,  A ,  involves  an  integral  over  a  rapidly 
oscillating  function  with  scale  length  e.  The  integral  over  such  a  function 
is  asymptotically  small  except  in  the  region  where  the  oscillations  become 
stationary.  Therefore,  coupling  between  the  modes  is  negligible  except  at 
the  stationary  point  where  the  integrand  in  the  exponential  vanishes.  The 
function,  A,  has  a  a  stationary  point  at  the  crossover  frequency,  u>cr,  where 
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d  =  0.  Because  the  contribution  to  the  overall  integral  in  (1.14)  is  exponen¬ 
tially  small  except  near  the  stationary  point,  we  may  expand  the  integral 
in  the  integrand  of  (1.14)  about  the  stationary  point  and  evaluate  it  using 
the  standard  methods  of  stationary  phase  [Bender  and  Orszag,  1978].  The 
overall  integration  may  then  be  extended  to  +oo  which  introduces  only  ex¬ 
ponentially  small  errors.  Then  we  find  that  the  coupling  coefficient,  up  to  an 
unimportant  phase,  is  given  by 


A(— oo)  = 


4e[r(a;cr)]1/2 


/+oo 

dsexp 

-OO 


/  -idfju^s2  \ 
\2t[r{ucr)Y12) 


(1.15) 


where 


df(u  —  ujcr) 


2<A^rJlclf2c2  (ui  +  ^2)2 

l,  -  i‘  n,n. 


(1.16) 


The  integral  may  then  be  rotated  in  the  complex  plane  in  which  case  it 
becomes  a  trivial  Gaussian  integral.  We  have  thus  obtained,  to  first  order  in 
k  the  asymptotic  form  of  the  coupling  to  the  LHCP  mode. 

By  comparing  the  Poynting  flux  along  the  magnetic  field 


z-S~z-E  xB~S(E*+E'+  + E*_E'_)  (1.17) 


of  the  incident  wave  above  the  coupling  region  to  that  of  the  coupled  wave 
below  the  coupling  region  we  may  obtain  the  coupling  coefficient 


C  =  b 


Z  *  S below 

7 r  k4 

Z  *  Sabot,e 

8e  yjrd' 

kA 

^Lek  II 

K\\ 


(1.18) 


This  result  may  be  naturally  extended  to  larger  values  of  k  if  we  consider 
the  embedded  dispersion  relation  discussed  in  Eq.  (1.7).  The  coupling  which 
occurs  at  the  crossover  frequency  arises  from  the  branch  point  which  is  the 
solution  to  Eq.  (1.11).  As  discussed  previously,  two  of  the  solutions  to  this 
equation  correspond  to  two  complex  conjugate  roots  which  occur  near  the 
crossover  frequency.  Moreover,  the  coupling  may  be  described  in  terms  of 
the  embedded  Eq.  (1.6)  with 


Id2  +  k4/A 

Q{z)  =  217^- 


(1.19) 
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The  function  d  has  a  zero  at  the  crossover  frequency,  so  that  locally  d  = 
dCr  +  d'cri f  with  dcr  =  0.  The  coupling  is  described  by  the  potential  equation 


2  2sr 


K 


(1.20) 


where  the  sum  function  has  been  evaluated  at  the  crossover  frequency  (the 
shape  of  the  potential  is  not  affected  too  much  by  the  expansion  of  s,  and  in 
particular  the  complex  turning  points  are  independent  of  the  form  of  s).  Eq. 
(1.20)  then  takes  the  form  of  a  Weber  equation.  The  transmission  coefficient 
for  this  equation  is  well  known  [Heading,  1962;  Lashmore- Davies  et  al.,  1985; 
Fuchs  et  al.,  1985]  and  takes  the  form 


T  =  exp 


(1.21) 


For  small  values  of  k,  T  «  1  -  C  as  obtained  in  Eq.  (1.18). 

As  is  clear  from  the  preceding  analysis,  the  coupling  occurs  as  the  result 
of  finite  k  and  increases  rapidly  for  small  angles  of  propagation.  For  larger 
k,  the  coupling  approaches  unity  which  means  that  a  complete  transfer  of 
power  from  the  predominately  RHCP  to  the  predominately  LHCP  mode 
occurs.  For  infinite  magnetic  field  scale  length  e  — y  0  so  that  coupling  is 
complete  except  at  k  =  0.  In  other  words,  waves  propagate  as  if  they  remain 
on  their  appropriate  dispersion  surface  and  do  not  couple  to  the  other  surface 
except  for  parallel  propagation.  In  that  sense,  the  scale  length  determines 
the  strength  of  the  coupling,  and  provides  a  criterion  for  the  validity  of  ray¬ 
tracing  in  the  plasma. 


3.2  Tunneling  Phenomena 

The  dispersion  relations  for  the  propagating  modes  of  the  differential  Eq. 
(1.1)  are  characterized  by  various  cutoffs  (k  -y  0)  and  in  particular  by  a 
cutoff-resonance  gap.  For  parallel  propagation,  the  modes  are  uncoupled 
and  the  LHCP  ion-cyclotron  exhibits  the  characteristics  of  the  “gap”  while 
the  RHCP  magnetosonic  branch  remains  essentially  constant.  The  LHCP 
branch  is  characterized  by  the  dispersion  relation  n_  —  l  where  /  has  a  zero 
at  the  cutoff  frequency,  ujco,  and  a  resonance  at  the  heavy  ion  gyrofrequency, 
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ficl.  For  non-parallel  propagation  the  situation  is  more  complex  in  that  the 
refraction  indices,  n±,  suffer  a  discontinuity  at  the  gyrofrequency  where  / 
diverges.  Across  this  discontinuity,  /  varies  as  1/u;  —  f)ci  which  changes  sign. 
Hence,  although  the  magnetosonic  solution  appears  to  be  continuous  across 
the  resonance,  there  is  in  fact  a  change  in  branch  n±  — >  The  existence  of 
such  a  discontinuity  signals  one  of  the  most  critical  aspects  of  the  problem 
which  we  are  considering  to  which  we  shall  return  many  times.  Namely,  the 
solution  is  not  continuous,  is  multivalued,  and  is  characterized  by  a  branch 
cut.  We  must  work  to  carefully  understand  both  the  physics  and  mathemat¬ 
ics  involved  in  the  multivalued  nature  of  the  solution  and  realize  that  the 
proper  choice  of  the  branch  of  this  solution  is  essential  to  our  analysis.  The 
proper  continuation  of  the  dispersion  relation  is  such  that  the  magnetosonic 
branch  is  not  continuous  across  the  resonance.  Hence,  although  the  ion- 
cyclotron  branch  apparently  mode  converts  into  the  magnetosonic  branch, 
the  asymptotic  behavior  of  the  ion-cyclotron  branch  below  the  resonance  is 
determined  primarily  by  the  asymptotic  behavior  of  the  ion-cyclotron  branch 
above  the  resonance. 

The  existence  of  a  “gap”  in  the  dispersion  relation  indicates  that  the  WKB 
solutions  for  a  particular  mode  are  cut  off  and  are  decaying  exponentials 
over  some  region  of  space.  In  general,  the  solutions  of  the  dispersion  relation 
across  this  “gap”  decay  such  that  only  exponentially  small  amplitude  may 
be  found  on  the  other  side  of  the  “gap.”  In  such  a  case,  the  mode  will  have 
been  reflected.  However,  if  the  “gap”  is  sufficiently  small  such  that  it  is  only 
a  fraction  of  a  wavelength  in  size,  substantial  tunneling  may  occur  in  which 
case  large  wave  amplitude  will  be  found  on  the  other  side  [ Budden ,  1985]. 

The  existence  of  a  pole  in  the  differential  equation  is  troublesome  in  that 
it  is  not  physically  correct,  and,  in  addition,  it  involves  the  complications 
of  a  regular  singular  point  which  generally  leads  to  multivalued  solutions  of 
a  differential  equation  except  under  extraordinary  circumstances.  As  far  as 
physical  correctness  is  concerned,  it  is  obvious  that  a  higher  order  expansion 
of  the  Z-functions  in  terms  of  the  thermal  velocity  will  introduce  corrections 
to  the  dispersion  relation  that  will  remove  the  singularity  (the  Z-functions 
are  well  defined  at  the  plasma  resonance,  hence  the  corresponding  differential 
equation  is  also  well  defined).  These  corrections  arise  in  the  form  of  a  higher 
order  terms  in  the  differential  equation  which  contain  the  essentials  of  the 
physics.  Physically,  the  pole  corresponds  to  dissipation  which  may  be  un¬ 
derstood  in  terms  of  strong  ion  cyclotron  damping.  In  order  to  understand 
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the  absorption,  we  may  consider  a  conserved  quantity  associated  with  the 
differential  equation  which  corresponds  to  the  Poynting  flux. 

From  Eq..  (1.1)  and  its  adjoint,  we  can  form  the  conserved  quantity 

J  =  -  tfttf')  (1.22) 

which  is  conserved  in  space  according  to 

J‘  =  2^f(M  -  Mf)ty  (1.23) 

where  denotes  real  derivatives  along  a  contour  with  fixed  imaginary  part 
(recall  that  dip* jdz  ±  (dip/dz)*  except  along  a  line,  z  —  x  +  iy,  with  y 
fixed).  The  conserved  quantity,  J,  corresponds  to  the  Poynting  flux  along 
the  2-axis,  z  *  (E  x  B),  and  is  conserved  along  the  real  axis  because  the 

functions  r  and  /  are  real.  The  function  /,  however  is  singular  at  the  heavy 

ion  gyrofrequency  so  that  we  expect  a  discontinuous  jump  in  the  value  of  J 
across  the  singularity.  The  simplest  model  which  approximates  the  functions 
r  and  /  near  the  resonance  and  retains  the  essential  characteristics  of  the 
dispersion  relation  is  one  in  which 

r  «  fcr2,  /  «  kf(  1  -a/z)  (1.24) 

where  a  linearized  form  of  the  magnetic  field,  B  —  £?q(1  —  2),  has  been  used. 
The  linearized  function,  /,  retains  the  resonance  and  cutoff,  while  variations 
in  the  function,  r,  are  not  substantial.  The  three  parameters  fcr,  k /,  and  a 
may  be  chosen  to  retain  the  location  of  the  crossover,  cutoff,  and  ion-ion 
frequencies  so  that  the  essential  characteristics  of  the  dispersion  relations  are 
obtained  with  quantitative  accuracy.  Near  the  resonance,  the  behavior  of  $ 
is  completely  dominated  by  the  singularity  in  /  so  that  higher  order  terms 
in  z  are  unimportant  for  determining  the  jump  in  J  across  the  cyclotron 
resonance. 

Near  the  gyrofrequency  (2=0),  the  cold  plasma  dispersion  relations  ex¬ 
hibit  an  unphysical  resonance.  If  we  were  to  incorporate  more  physical  con¬ 
siderations  into  our  analysis,  such  as  the  thermal  effects,  we  would  find  that 
the  solution  to  the  dispersion  relation  near  the  gyrofrequencies  is  of  the  form 
co  ~  f]cj  —27,7  >  0,  such  that  the  mode  is  strongly  damped.  In  light  of 
our  approximation  for  the  magnetic  field  the  actual  pole  should  be  at  a  lo¬ 
cation  u  —  ( flci  —  ij)  =  u  —  u(  1  —  z)  27  =  0,  that  is,  at  the  location 
*  =  —iy/u  =  -iy  which  lies  below  the  real  axis  rather  than  at  2  =  0. 
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We  may  determine  the  discontinuity  across  the  resonance  by  considering 
the  change  in  J  along  the  path  2  =  x±iy.  From  Eq.  (1.23),  the  discontinuity 
is  given  by 

Ay  =  ±2kfa  f  2  V  (1.25) 

y  1  J- 00  x2  +  y 2 

where  (j)  the  lower  component  of  For  small  2,  <f>  behaves  like  the  Hankel 
functions  or  integrals  thereof,  which  are  regular  [Abramowitz  and  Stegun, 
1964],  so  that  the  discontinuity  is 

lim  As  =  ^2nkfa\(j)\l=0  =  T A,  A  >  0  (1-26) 

y-+0 

and  represents  energy  which  is  absorbed.  That  is  <  J  where  the 
superscript  refers  to  the  value  along  the  positive  (negative)  real  axis.  If  there 
is  no  flux  in  or  out  on  the  positive  real  axis  (J+  =  0)  there  is  a  positive 
energy  flow  from  the  left  (<!T~  >  0),  which  is  absorbed.  On  the  other  hand, 
if  there  is  no  flux  from  the  left  (J~  =  0),  then  there  is  a  negative  flux  from 
the  right  (J+  <  0)  sending  wave  energy,  which  is  absorbed,  into  the  system. 
The  value  of  this  energy  loss  clearly  depends  on  the  details  of  the  boundary 
conditions  in  that  the  solution  xp  must  be  known  in  order  to  calculate  the 
energy  loss.  Physically,  we  may  interpret  this  energy  loss  as  resulting  from 
resonant  absorption  at  the  gyrofrequency.  Such  absorption  typically  occurs 
because  the  cold  plasma  wave  mode  converts  into  a  thermal  wave  which  is 
readily  absorbed.  Conversely,  if  the  pole  lies  above  the  real  axis,  a  similar 
procedure  results  in  an  energy  gain  which  is  not  physical. 

We  may  obtain  a  substantial  amount  of  information  concerning  the  asymp¬ 
totic  behavior  of  the  solutions.  The  WI<B  solutions  in  this  case  are  of  the 
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and  the  zero  order  sum  and  difference  functions  are 


*  =  *+*  >0 

(1.30) 

so  that  the  four  WKB  solutions  take  the  form 

0+  =  exp(±i/e(n+z  —  /?+  log  z/2)) 

Qt  =  exp(±i/e(n„z  -  /?_  log  z/2)) 

(1.31) 

where 

k?a±$  +  \fS2  +  k4/ 4 

P±  =  - r . . 

2  n±\A2  +  K4  j\ 

(1.32) 

and  the  coupling  is  characterized  by  the  multivalued  nature  of  the  logarithm. 

From  these  expressions  we  can  determine  the  transmission  coefficients  for 
the  case  of  incident  magnetosonic  waves  or  ion-cyclotron  waves.  We  may 
extend  the  WKB  solutions  into  the  complex  plane  and  use  the  dominance  of 
the  various  solutions  to  determine  transmission  properties  [ Heading ,  1962]. 
The  solutions  on  the  positive  real  axis  are  of  the  form 

V’  =  A0+  +  BQ_  -f  CQt  +  F*0+,  Arg(z)  =  0  (1.33) 

where  we  have  written  them  in  order  of  dominance  in  the  upper  half  complex 
plane.  The  lower  component,  < j>,  of  is  of  the  same  form  with  coefficients 
which  are  related  to  A,  B,  C,  and  D  by  the  dispersion  relation.  The  WKB 
approximations  in  the  complex  plane,  however,  are  at  best  limited  to  validity 
within  certain  sectors  due  to  the  Stokes  phenomenon  which  results  from  the 
fact  that  the  exponential  solutions  are  too  simple  to  adequately  describe  the 
qualitative  behavior  of  the  actual  functions.  Although  the  coefficient  of  the 
dominant  solution  within  a  given  region  remains  unchanged,  the  subdomi¬ 
nant  solutions  are  liable  to  change.  Moreover,  on  the  grounds  of  continuity 
we  expect  that  if  the  coefficient  of  the  more  dominant  solutions  vanish,  the 
coefficient  of  the  most  dominant  remaining  solution  is  unchanged.  To  accom¬ 
modate  this  change,  we  make  a  discontinuous  adjustment  of  the  subdomi¬ 
nant  terms  along  the  Stokes  lines  where  the  subdominant  terms  are  least 
significant.  The  adjustment  should  be  proportional  to  the  amplitude  of  the 
dominant  solution  so  that  it  will  vanish  if  the  dominant  solution  vanishes. 
The  Stokes  lines  in  this  case  are  along  the  positive  imaginary  axis  while  the 
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anti-Stokes  lines  are  along  the  positive  real  axis.  According  to  the  Stokes 
phenomenon  the  solution  takes  the  general  form 

0=  AQ+  +  (B-AS1)QZ  +  {C-AS2-BS3)Ot 

+(D  -  AS4  -  BS5  -  CS6)0+,  Arg(z)  =  -n  (1.34) 

on  the  negative  real  axis.  We  have  added  Stokes  constants  across  the  region 
of  least  dominance  according  to  a  reasonable  prescription  which  retains  the 
continuity  of  the  equations  in  the  absence  of  the  more  dominant  solutions. 
The  boundary  conditions  for  downgoing  ion-cyclotron  waves  are  given  by 

>4  =  1,  B  =  Si,  C  =  D  =  0  (1.35) 


while  for  downgoing  magnetosonic  waves  the  boundary  conditions  are 


>4  =  0,  5  =  1,  C  =  D  =  0 


(1.36) 


We  may  immediately  write  down  the  asymptotic  form  of  the  transmission  co¬ 
efficients  in  terms  of  the  logarithmic  singularity.  The  leading  order  coefficient 
for  downgoing  ion-cyclotron  waves  is  given  by 


ry,  l®+l>lrg(z)=7r _ -r0+/t 

T,c  =  ]eji1  -  e 


(1.37) 


\Arg(z)=0 


whereas  for  downgoing  magnetosonic  waves  in  the  absence  of  an  incident 
ion-cyclotron  wave  has  no  transmitted  flux  in  the  ion-cyclotron  mode  and 
transmission  to  the  magnetosonic  mode  given  by 


Tms  = 


1^—  l>4r-ff(z)=7r  _  e_7r 0-/c 
l®-l^rj(z)=0 


(1.38) 


We  can  make  a  number  of  observations  concerning  these  two  results. 
First,  for  k  «  d,  we  find  from  Eq.  (1.32)  that  Tic  is  simply  the  Budden 
coefficient  described  previously.  That  is  (3+  =  fc/a(  1  —  0(k4))  which  has  little 
dependence  on  k,  so  that  transmission  for  this  mode  may  be  considered  to 
be  essentially  Budden  tunneling  and  should  have  that  character.  With  in¬ 
creasing  k,  the  tunneling  diminishes.  For  large  /c,  we  find  that  the  coefficient 
approaches  another  constant,  namely 

*  -  Wo  (1'39) 
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indicative  of  Budden  tunneling  for  nearly  perpendicular  propagation  near  the 
ion-ion  hybrid  frequency. 

On  the  other  hand,  for  the  downgoing  magnetosonic  wave  we  find  that 
for  small  k  Eq.  (1.32)  gives  an  intriguing  result.  The  transmission  coefficient 
then  takes  the  form 

Tms  =  <=-*-'«  *  exp  (1.40) 

We  may  relate  this  transmission  coefficient  to  the  coupling  coefficient  ob¬ 
tained  from  the  saddle  point  analysis.  If  we  recognize  that  difference  function 
is  of  the  form 

r  kfa 

d=-s+~t  <L4i) 

we  find  that 

2P 

dcr  =  Jfa  (L42) 

which  then  demonstrates  that  the  transmission  coefficient  is  the  same  as 
the  coupling  coefficient  established  in  Eq.  (1.18).  We  may  make  a  physical 
interpretation  of  this  result.  For  parallel  propagation  transmission  of  the 
magnetosonic  mode  is  complete  (T  =1).  However,  coupling  increases  rapidly 
with  larger  values  of  /c,  and  we  find  that  the  transmission  in  the  magnetosonic 
mode  falls  off  substantially. 

The  values  of  the  Stokes  constants  may  not  be  unambiguously  determined 
using  the  methods  of  the  previous  section  although  a  number  of  relationships 
between  the  amplitudes  of  the  various  waves  may  be  determined.  For  the 
case  of  a  “large”  gap,  however,  we  may  actually  determine  the  coefficients 
unambiguously.  In  this  case  the  WKB  solution,  0^,  has  negligible  amplitude 
at  —oo,  i.e.  /?+  ->  o o  so  that  there  is  no  transmission  in  that  mode.  Fur¬ 
thermore,  the  restriction  that  no  upgoing  ion-cyclotron  wave,  0+,  be  found 
at  — oo  implies  that  A—  >  0  so  that  no  absorption  occurs.  In  this  case,  after 
substantial  and  lengthy  algebra  detailed  in  Appendix  A,  it  is  possible  to  show 
that  in  the  limit  of  a  large  “gap”  the  transmission  and  reflection  coefficients, 
as  defined  in  the  appendix  to  be  the  ratio  of  the  Poynting  flux  of  a  particular 
wave  to  that  of  the  incident  wave,  are  given  by 

Cjc  =  1  -  =  1  -  C 

Rlc  =  e-2^-A  =  C2  (1.43) 

Cfc  =  e“^-A(l  -  e-^-/€)  =  C(1  -  C) 
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Figure  3:  Illustration  of  the  meaning  of  the  coefficients  for  a  downgoing  ion 
cyclotron  wave.  Waves  are  reflected  back  along  the  two  upward  propagating 
modes  while  a  wave  is  transmitted  on  the  downgoing  magnetosonic  branch. 

where  the  meaning  of  the  coefficients  are  illustrated  in  Figure  3.  For  a 
magnetosonic  wave  incident  from  above  we  find  that 

Tms  =  e-^-/£  =  C 

Rms  =  (1  -  e-^-/')2  =  (1  -  C)2  (1.44) 

C«s  =  e-^-/£(l  -  e-^-R)  =  C{\  -  C) 

This  process  is  again  illustrated  in  Figure  4. 

As  stated  previously,  in  the  case  of  large  gap,  wave  propagation  is  entirely 
determined  by  the  interaction  at  the  crossover  frequency  which,  as  we  have 
just  shown,  is  characterized  by  the  coefficient  C  =  e~nl3~R.  In  order  to 
understand  these  results  in  an  intuitive  manner,  consider  the  situation  of 
Figure  3.  A  downgoing  ion-cyclotron  wave  with  unit  flux  encounters  the 
crossover  frequency  at  which  point  a  fraction,  C,  is  transferred  from  branch 
II  to  branch  I.  A  fraction,  1  —  C,  continues  down  along  mode  II  and  is  the 
coupled  transmission  coefficient  Cjc.  The  wave  which  couples  to  branch 
I  reflects  at  the  cutoff  frequency  and  once  again  encounters  the  crossover 
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Figure  4:  Illustration  of  the  meaning  of  the  coefficients  for  a  downgoing 
magnetosonic  wave.  Waves  are  reflected  back  along  the  two  upward  prop¬ 
agating  modes  while  a  wave  is  transmitted  on  the  downgoing  magnetosonic 
branch. 

frequency  as  an  upgoing  wave.  This  reflected  wave  has  amplitude  C.  Upon 
coupling,  a  fraction,  C,  of  this  reflected  wave  flux  is  coupled  to  branch  II 
comprising  the  reflection  coefficient  RIC  =  C2.  The  fraction,  1  -  C,  remains 
on  branch  I  and  is  the  coupled  reflected  wave  Cfc  —  C(  1  —  C).  Although 
in  the  limit  of  a  small  “gap”  the  transmission  coefficient  Tjc  is  a  Budden 
tunneling  term,  in  the  limit  of  a  large  “gap”  TIC  vanishes. 

The  case  of  an  incident  magnetosonic  wave  is  illustrated  in  Figure  4. 
In  this  case,  a  magnetosonic  wave  of  unit  amplitude  is  incident  along  branch 
I.  Coupling  occurs  at  the  crossover  frequency  and  a  fraction,  C,  of  the  inci¬ 
dent  wave  is  coupled  to  branch  II  and  comprises  the  transmission  coefficient, 
Tms  —  C .  The  fraction,  1  —  C,  remains  on  branch  I,  is  reflected  at  the  cutoff 
frequency,  and  again  encounters  the  crossover  frequency  as  an  upgoing  wave. 
The  amplitude  of  this  reflected  wave  is  1  —  C .  Again  coupling  occurs  at 
the  crossover  frequency  and  a  fraction,  C,  of  the  reflected  wave  is  coupled 
onto  mode  II  giving  the  coupled  reflected  wave  C%JS  =  C(1  -  C ),  and  the 
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remaining  wave  power,  Rms  =  (1  —  C )2  is  reflected  along  branch  I.  There  is 
no  coupled  mode  CjjS  in  the  limit  of  large  “gap.”  As  one  can  see,  for  the 
case  of  a  large  frequency  gap,  the  coefficients  are  well  described  in  terms  of 
the  coupling  at  the  crossover  frequency. 

In  the  case  of  a  small  “gap”  it  is  not  possible  to  determine  the  Stokes 
parameters  unambiguously  although  a  number  of  relationships  between  them 
may  be  derived  using  symmetry  properties.  In  the  following  chapter  we  will 
determine  the  coefficients  numerically.  The  results  which  we  have  obtained, 
however,  are  useful  in  that  they  provide  both  physical  insight  and  a  check  of 
our  numerical  procedure.  The  method  which  we  have  used  for  the  Budden 
equation  is,  however,  sufficient  to  determine  the  coefficients  for  large  k  in 
which  case  only  one  mode  propagates  at  ±oo. 

An  alternate  approach  which  makes  use  of  the  Fourier  transform  of  Eq. 
(1.1)  has  been  presented  by  Le  Queau  and  Roux  [1992].  Their  approach  is 
based  upon  the  asymptotic  expansion  of  the  integral  transform  of  the  ap¬ 
proximate  solution  to  the  differential  equation.  The  solution  is  continued 
into  the  complex  plane  much  in  the  manner  of  the  phase  integral  approach, 
and  they  derive  the  same  reflection,  coupling  and  transmission  coefficients 
as  detailed  above.  Their  result  expands  upon  the  approach  detailed  above 
and  is  useful  in  the  regime  in  which  the  “gap”  is  small.  In  the  limit  of  a 
large  “gap,”  the  results  obtained  above  are  in  agreement  with  their  calcula¬ 
tions.  Furthermore,  their  results  are  consistent  with  the  relations  that  may 
be  derived  between  the  coefficients.  Our  work  substantiates  the  coefficients 
which  they  obtained  in  an  independent  manner.  The  reader  is  referred  to 
their  paper  for  the  complete  set  of  coefficients. 

3.3Phase  Integral  Solution  for  Large  k — Two  Propa¬ 
gating  Modes 

As  one  can  see  from  Eq.  (1.28),  for  large  enough  a  <  k2,  so  that  two 
of  the  solutions  become  complex.  That  is,  we  find  that  n+  is  real  and  n_ 
is  imaginary.  This  means  that  two  of  the  asymptotic  solutions  are  growing 
along  the  real  axis  (a  notably  unphysical  situation).  When  this  is  the  case, 
the  positive  and  negative  real  axes  become  Stokes  lines  for  two  of  the  solutions 
while  the  positive  and  negative  imaginary  axes  remain  the  Stokes  lines  for 
the  other  two  solutions.  The  anti-Stokes  lines  are  to  be  found  asymptotically 
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along  the  rays 

n“r  =  ±n°_y  (1.45) 

where  2  =  x+ iy  is  the  complex  coordinate.  This  situation  is  easier  to  evaluate 
because  each  solution  is  dominant  in  one  of  four  sectors  in  the  complex  plane. 

At  this  point,  we  need  to  consider  the  concept  of  Stokes  phenomena  in 
somewhat  more  detail.  Tracing  Stokes  coefficients  in  reality  is  only  valid  in 
the  sense  that  in  the  region  of  maximum  dominance,  the  coefficient  of  the 
exponential  term  must  agree  with  the  exact  solution  in  an  asymptotic  sense. 
However,  we  actually  require  further  information  concerning  the  asymptotic 
behavior  of  the  subdominant  terms  in  order  to  extract  the  reflection  and 
transmission  coefficients.  This  information  is  important  because  we  must 
determine  the  value  of  the  oscillatory  solutions  along  the  real  axis  in  the 
presence  of  an  unphysical  growing  mode. 

In  previous  work  [  Johnson ,  1992]  we  demonstrated  using  an  integral  equa¬ 
tion,  that  the  Stokes  coefficients  for  the  asymptotic  solutions  of  the  fourth 
order  equation  which  we  are  considering  change  abruptly  at  definite  values 
of  Arg(z)  in  the  region  of  maximum  subdominance.  It  would  appear  that 
these  jumps  occur  at  three  distinct  values  of  Arg(z) — along  the  Stokes  lines 
and  the  bounding  anti-Stokes  lines  of  the  maximally  subdominant  region. 
The  change  only  occurs  in  the  coefficient  of  the  most  subdominant  solution. 
We  may,  therefore,  determine  the  appropriate  asymptotic  behavior  of  the 
solution  of  interest  if  the  relevant  WKB  solutions  involved  in  the  analysis 
are  not  maximally  subdominant  in  the  region  where  they  must  be  evaluated. 
Because  the  solutions  with  growing  (decaying)  exponential  behavior  do  not 
contribute  to  the  Poynting  flux  on  the  real  axis,  we  need  only  evaluate  the 
coefficients  of  the  oscillatory  solutions  in  order  to  determine  the  transmission 
and  reflection  coefficients.  The  oscillatory  solutions  are  not  maximally  sub¬ 
dominant  along  the  real  axis,  and  therefore,  we  may  determine  the  coefficients 
unambiguously.  In  order  to  incorporate  the  Stokes  phenomena  we  introduce 
the  Stokes  multipliers  such  that  upon  crossing  a  Stokes  line  on  which  a  WKB 
solution  is  maximally  subdominant  the  coefficient  of  that  term  equals  the  co¬ 
efficient  before  crossing  the  Stokes  lines  +  J2j  SjX  the  coefficients  of  the  jth 
dominant  term  on  the  Stokes  line  [ Heading ,  1957;  Buddert,  1985]. 

Appropriate  boundary  conditions  require  that  at  ±00  the  exponentially 
growing  modes  be  absent.  In  addition,  appropriate  boundary  conditions 
require  that  no  upgoing  wave  be  found  at  —00.  In  order  to  satisfy  these 
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Figure  5:  Sectors  of  the  complex  plane  divided  according  to  regions  of 
dominance  and  subdominance  of  the  asymptotic  WKB  solutions.  The  Stokes 
lines  and  anti-Stokes  lines  are  solid  and  dashed  respectively. 


conditions,  we  assume  along  or  above  the  positive  real  axis  that  the  solution 
takes  the  form 


ip  =  AQ+  +  BQ_  +  CQt.  +  Z)0+,  Region  II 

(1.46) 

Using  the  rules  stated  above  we  may  continue  the  solution  into  the  following 
sectors  illustrated  in  Figure  5. 

rP  =  AQl  +  BQZ  +  C’Qt  +  DQ+,  Region  III 

(1.47) 

with 

C’=C-ASi-  BS2  -  DS4 

(1.48) 

4’  =  A’Q-  +  BOZ  A  C'et  +  Del,  Region  IV 

(1.49) 

with 

A'  =  A-  BU2-  CU3  -  DU4 

(1.50) 
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and 

ip  =  A0+  +  BQ_  +  cet  +  Z/0+,  Region  I  (1.51) 

with 

D*  =  D  +  AV1  +  BV2  +  CV3  (1.52) 

where  S,  U,  and  V  are  Stokes  constants.  In  order  to  satisfy  the  boundary 
conditions  for  downgoing  waves,  we  require  that 


B  =  C=  D'  =  0,  A  =  1 


(1.53) 


Then  we  may 
ficients 


immediately  write  down  the 


TIC  = 


|©-12 


+  IArg(z)=pi 


|©-'2 


+  IArg(z)=0 


transmission  and 
=  e_7r/3+/£ 


reflection  coef- 
( 1 .54) 


and 


Ric 


Urg(e)=0  |T/  |2 
Terp  l^il 

+  \Arg{z)=0 


(1.55) 


In  order  to  determine  the  Stokes  constant  V\  we  may  consider  the  situation 
with  B  =  C  =  D  =  0,A  =  l  from  which  we  can  compute  the  Stokes  constant 
if  we  again  invoke  the  conservation  law.  As  it  turns  out,  on  the  real  axis, 
the  exponentially  growing  solutions  do  not  contribute  to  the  Poynting  flux 
so  that  the  conserved  quantity  J/n+  is  found  to  be 


1  Arg(z)  =  0 

e^+A  Arg[z )  = 

e-^+A-lVilV^/*  Arg{z)  =  v 


Once  again,  from  the  energy  discontinuity 


-A  =  1  -  e^+A  =  -1  +  -  J Vi |2eTr/?+  /c 

we  obtain  the  Budden-like  reflection  coefficient 


Ric  =  l^il2  =  (1  -  T/C)2 


(1.56) 


(1.57) 


(1.58) 
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3.4  Mode  Coupling  at  Low  Frequencies,  Large  k 


The  result  which  we  have  just  obtained  for  large  values  of  k  does  not  depend 
critically  on  the  crossover  frequency  in  that  the  factor  (3+  depends  weakly 
on  the  location  of  the  crossover  frequency.  From  the  analysis  described  in 
Section  3.1,  we  find  that  for  large  values  of  «,  the  saddle  point  associated 
with  the  crossover  frequency  is  located  far  from  the  real  axis  and  is  near  the 
hydrogen  cyclotron  frequency.  As  a  result,  coupling  is  not  strongly  dependent 
upon  the  location  of  the  saddle  point  associated  with  the  crossover  frequency. 
However,  in  Figure  2,  we  find  that  the  coupling  point  associated  with  zero 
frequency  at  small  k,  approaches  the  oxygen  resonance  for  large  values  of  k 
and  is  therefore  of  primary  interest.  As  described  previously,  the  transition 
from  four  propagating  modes  to  two  propagating  modes  occurs  at  a  critical 
value  of  k.  Assuming  that  r  and  /  take  a  simple  form  such  as  in  Eq.  (1.24), 
the  cutoff  condition  (1.9)  implies  that 


1  -  «2/4fcr2 

co  °i  -<rKyk*kf  (L59) 

We  then  find  that  the  abrupt  change  in  the  behavior  of  the  modes  occurs  at 
the  critical  value 

A-2i-2 

(1.60) 


h2h2 

rC  — 


For  values  smaller  than  this  value  of  k,  the  mode  conversion  processes  are 
determined  primarily  from  the  coupling  point  associated  with  the  crossover 
frequency.  For  intermediate  values  of  «,  k2kf/cr  <  k2  <  2 kr,  four  modes 
propagate  in  the  vicinity  of  the  heavy  ion  resonance  although  the  magne- 
tosonic  branch  is  cut  off  at  some  point  below  the  gyrofrequency  and  does  not 
propagate  above  the  gyrofrequency.  In  this  case,  the  dominant  saddle  point 
is  that  associated  with  zero  frequency  at  small  k  and  the  oxygen  resonance 
at  large  k.  For  larger  k,  k2  >  2 kr,  the  downgoing  ion-cyclotron  wave  is  cutoff 
above  the  oxygen  resonance  at  a  frequency  which  approaches  the  ion-ion  hy¬ 
brid  frequency  and  only  two  modes  propagate.  In  this  case  the  transmission 
properties  may  be  obtained  by  expanding  about  the  coupling  point  associ¬ 
ated  with  the  cutoff  condition  rl  —  k2s  =  0.  In  the  ensuing  analysis  we  will 
follow  the  reasoning  of  Lashmore-Davies  et  al.  [1985]  and  Fuchs  et  al.  [1985], 
and  we  will  consider  the  cases  of  intermediate  and  large  k  separately. 

The  expression  for  the  location  of  the  saddle  point  in  the  complex  plane 
is,  as  before,  given  by  Eq.  (1.11)  As  discussed  in  Section  3.1,  the  solutions  to 
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this  condition  are  complex  conjugate  pairs  which  satisfy  a  quartic  equation 
as  plotted  in  Figure  2.  In  this  case  the  potential  again  takes  the  form 
(1.19).  The  fact  that  the  roots  occur  in  complex  conjugate  pairs  means  that 
the  function,  d 2  +  k4/4,  will  have  the  form  of  a  complex  conjugate  barrier 
equation  described  by  the  Weber  equation.  It  is  once  again  both  tractable 
and  somewhat  reasonable  to  assume  the  form  of  r  and  /  as  assumed  in  Eq. 
(1.24).  In  this  case,  the  potential  takes  the  form 


0(0 


S2  +  *74 
4(2(s  -  *72) 


{((  -  7)2  +  (?) 


(1.61) 


where 


_  kfa  S  ±  in2  / 2 

=  ~2~£2  +  /c74 


(1.62) 


is  the  location  of  the  complex  conjugate  roots  (for  large  k  the  imaginary  part 
dominates  as  is  the  case  for  the  actual  solutions  to  the  third-order  equation). 
Once  again  the  essential  asymptotic  features  of  the  embedded  dispersion 
relation  are  retained  if  we  evaluate  the  constant  factor  in  Eq.  (1.61)  at  some 
convenient  point,  such  as  the  real  frequency  above.  If  we  then  make  the 
judicious  choice  for  the  variables 


j  2 _  82  +  k4/4 

~°  "  er(Sr  ~  K*/2) 


and 


kfa 


(3=-  =  - 

2  2(52  +  k4/4)  4 


(1.63) 


(1.64) 


(where  the  index  r  indicates  that  the  function  is  to  be  evaluated  at  £r)  then 
the  embedded  equation  is  of  the  form  of  the  Weber  equation  which  has  the 
transmission  coefficient  [, Lashmore- Davies  et  al. ,  1985;  Fuchs  et  al,  1985] 


J1  _  e-2irfc0j92A 

where  the  argument  of  the  exponential  takes  the  form 

k,a/45 


he2  =  7 


8  y/fS2  +  k4/4)(/c4/4£  +  «:2/2  -  kf) 


(1.65) 


(1.66) 
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Near  the  upper  limit  of  validity,  the 
limit.  For  k2/ 2  =  k 2 


whereas  the  equivalent  coefficient  from  the  saddle  point  becomes 


two  results  have  a  reasonably  simple 

*?“  n  e-n 


2  kaf  =  p^  (1.68) 

yielding  roughly  the  same  expression.  Over  the  range  of  densities  for  which 
absorption  is  strong  (0.5-5%),  the  ratio,  k2/2y/8  takes  values  ranging  from 
1.3  to  0.6  depending  on  the  density  ratio  and  the  two  results  are  in  reasonable 
agreement.  Hence,  we  can  understand  that  the  coupling  which  arises  at  low 
frequencies  results  from  coupling  at  the  saddle  point  which  occurs  near  the 
oxygen  resonance  for  large  values  of  k. 

The  form  of  the  reflection  coefficient  may  be  understood  physically  in 
the  same  manner  that  we  understood  coupling  for  smaller  k.  In  this  case, 
coupling  between  the  modes  occurs  below  the  heavy  ion  gyrofrequency.  As 
illustrated  in  Figure  6,  incident  waves  along  mode  II  couple  at  the  complex 
coupling  point.  A  fraction,  Tjc  —  T  =  e-,r^+ ,  is  transferred  to  branch  III.  The 
remaining  wave  energy,  1  —  T,  continues  along  branch  II,  is  reflected  at  the 
cutoff  frequency  for  mode  II  and  once  again  encounters  the  coupling  point. 
The  two  waves  couple  again  and  a  fraction,  1  —  T,  of  the  remaining  wave 
power  is  found  on  branch  II  giving  the  reflection  coefficient,  Ric  =  (1  —  T)2. 
The  rest  of  the  wave  power  is  absorbed  at  the  resonance. 

For  larger  values  of  k  we  must  expand  about  the  cutoff-resonance  pair 
which  arises  near  the  ion-ion  hybrid  frequency.  In  this  case  the  coupling 
potential  in  (1.6)  takes  the  form 


Q  = 


rl  —  k2s 
r  +  /  -  k2 


(1.69) 


Once  again,  the  functions  may  be  expanded  about  the  gyrofrequency,  fici  = 
cj(1  —  z).  After  suitable  definitions  and  algebraic  manipulation  detailed  in 
Appendix  B,  the  potential  is  found  to  take  the  Budden  form 


Q  =  ko(l  -  alz) 


(1.70) 
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Figure  6:  Illustration  of  the  physical  meaning  of  the  coefficients  for  a  down¬ 
going  ion-cyclotron  wave  for  large  angle  of  propagation.  Waves  are  transmit¬ 
ted  onto  mode  III  and  reflected  back  along  mode  II. 


The  transmission  coefficients  for  the  Budden  equation  are  well  known 

T  =  e~nkoQ/e,  R  =  (\-Tf  (1.71) 

In  the  limit  of  large  k  the  potential  takes  the  form 


Q  K  s  = 


*4  (^-f^X^-r^) 


(1.72) 


which  exhibits  a  cutoff  at  the  ion-ion  hybrid  frequency  and  a  resonance  at 
the  oxygen  gyrofrequency.  Expanding  s  about  the  gyrofrequency 


s  «  cr(l  —  kfa/2(rz) 

we  find  that  the  coefficient  takes  the  form 


k0a  =  (3+  = 


kfa 


(1.73) 


(1.74) 
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in  agreement  with  (1.39)  of  the  previous  section. 

Comparing  both  this  result  and  the  result  for  intermediate  we  find  that 
both  coefficients  are  of  the  Budden  form  with  transmission  and  reflection  co¬ 
efficients  given  by  T  and  (1  —  T )2  respectively.  These  results  compare  well 
with  the  form  (1.54,  1.58)  which  we  obtained  using  the  phase-integral  analy¬ 
sis.  In  this  analysis,  we  have  three  parameters:  kr,  ki  and  a.  These  three  pa¬ 
rameters  are  sufficient  to  retain  the  location  of  the  resonance  (flci),  crossover 
frequency  (o>cr),  cutoff  frequency  (u>c0),  and  ion-ion  hybrid  frequency  (u >,,). 
Thus,  in  some  sense  we  can  form  a  uniform  approximation  which  retains  the 
characteristics  of  the  crossover  and  cutoff  frequencies  for  small  k  but  retains 
the  characteristics  of  the  ion-ion  hybrid  resonance  for  perpendicular  propa¬ 
gation.  This  linear  approximation  is  good  if  the  coupling  is  dominated  by 
one  or  the  other  of  the  two  coupling  points,  but  if  both  coupling  points  are  of 
significant  importance,  one  suspects  that  the  linear  approximation  will  fail 
for  the  phase-integral  approach. 


4  Summary 

In  this  report  we  have  explored  the  analytic  properties  of  the  mode  con¬ 
version  processes.  We  have  found  that  near  the  crossover-cutoff-resonance 
triplet  the  propagating  modes  couple  strongly,  and  substantial  energy  may 
be  exchanged  between  modes.  For  small  values  of  k,  we  have  determined 
that  the  ion-cyclotron  mode  tunnels  in  a  manner  described  by  the  Budden 
equation,  and  substantial  wave  power  will  be  transmitted  and  absorbed  for 
nearly  parallel  propagation.  We  have  also  shown  that  this  tunneling  factor 
drops  off  with  larger  /c,  vanishing  as  k  becomes  the  order  of  1  (in  which  case 
the  approximation  is  not  valid).  For  larger  angles  of  propagation,  this  mode 
continues  to  be  characterized  by  Budden  tunneling;  however,  in  this  case, 
the  mode  does  not  interact  with  the  nonpropagating  magnetosonic  wave. 
For  large  angles  of  propagation,  the  mode  is  described  by  a  simple  expres¬ 
sion  for  tunneling  at  the  ion-ion  hybrid  resonance.  Because  the  “gap”  which 
characterizes  this  tunneling  is  smaller  than  that  for  parallel  propagation,  we 
expect  that  the  tunneling  and  absorption  will  be  of  the  same  order  as  for 
parallel  propagation. 

Incident  magnetosonic  waves  are  only  of  interest  for  reasonably  small  k 
because  for  larger  k  they  do  not  propagate.  For  parallel  propagation  the 
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coupling  is  negligible  and  the  mode  propagates  without  much  interaction. 
However,  the  coupling  increases  rapidly  as  k4.  For  large  k  the  coupling  is 
reasonably  complete  and  all  the  wave  power  is  reflected.  The  reflection  occurs 
because  the  magnetosonic  mode  couples  completely  to  the  LHCP  wave  which 
is  cut  off  far  above  the  oxygen  resonance  so  that  the  tunneling  barrier  becomes 
insuperable. 

In  the  rest  of  this  report,  we  explore  the  numerical  value  of  these  coef¬ 
ficients.  In  particular,  we  will  solve  the  differential  equations  numerically 
in  the  regimes  where  we  have  not  obtained  the  coefficients.  As  we  shall 
demonstrate,  the  basic  features  of  the  analytical  results  are  retained. 


5  Appendix  A:  Coefficients  for  a  Small  “Gap” 

In  this  appendix  we  consider  the  linear  approximation  in  the  limit  of  a  small 
“gap.”  Given  the  linear  approximations  (1.24),  the  solutions  of  (1.1)  may  be 
decomposed  into  the  appropriate  WKB  solutions  (1-27)  far  from  the  coupling 
region.  Along  the  positive  (negative)  real  axis,  the  solutions  are  of  the  form 
(1.33)  and  (1.34)  respectively.  A  large  energy  gap  implies  that  A  =  0  if  the 
ion-cyclotron  wave  is  not  found  below  the  coupling  region.  In  the  asymptotic 
regime  along  the  real  axis,  the  conserved  flux  is  of  the  form 


/+|A|2e-^  +/_|B|2e-’>-  — /_|C|V-  -f+\D\2e'+ 


while  along  the  positive  real  axis,  it  takes  the  form 


(1.A1) 


f+\A\2+f_\B+S1A\2-f_\C+S2A+S3B\2-f+\D+S4A+S5B+S6C\2  (1.A2) 


where  the  constants  f±  are  factors  which  relate  the  lower  and  upper  compo- 
nents,  xp  and  0,  of  the  differential  Eq.  (1.1) 


fdt  71  ± 


and 

V±  =  *P± 

A  large  “gap”  implies  that  7y+  is  very  large. 


(1.A3) 

(1.A4) 
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5.1  Magnetosonic  Wave 

The  boundary  conditions  for  a  downgoing  magnetosonic  wave  are 

,4  =  0,  £=  1,  (7  =  0,  D  =  0  (1.A5) 

for  which  we  find 

f  -  /-|S3|2  -/+|Ss|2  Arg(z)  =  0 
1  f_e  v~  Arg(z)  =  n 

Because  A  =  0  the  ratios  of  these  terms  correspond  to  the  coefficients.  Thus, 

Tms  =  e-"-,  Rms  =  |5s|J,  Cfcs  =  |t|S5|J  (1.A7) 

A  suitable  choice  of  three  other  boundary  conditions  is  sufficient  to  determine 
the  Stokes  coefficients,  IS3I,  |*?5|  and  |Se|.  We  find  that 

12  f- 


156,2  =  U{eV~  ~ 1} 
|S3|’ =  (e-'- -  If 


so  that 


Tms  =  e-*1-,  Rms  =  (1  -  e’"-)2,  C^s  =  e-”-(l  -  e"’1") 

5.2  Incident  Ion-Cyclotron  Wave 

The  boundary  conditions  for  an  ion-cyclotron  wave  are 

A  =  1,  B  =  Si,  (7  =  0,  D  =  0 

for  which  we  find 

f  f+-f_\s2  +  S',  S3|2  -  f+\s4  +  S,S5|2  Arg(z)  =  0 
J  \  /_|5i|2e",J-  Arg(z)  =  tt 


(1.A8) 


(1.A9) 


(1.A10) 


(1-All) 
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where  it  is  to  be  noted  that  we  have  ignored  the  term  with  the  factor  e-”+ 
along  the  negative  real  axis  in  that  it  is  exponentially  small.  We  may  then 
identify  the  coefficients  to  be  the  ratios  of  the  downgoing  fluxes  to  the  incident 
flux 

Cic  —  V~ ’  Cfc  =  Ric  =  \S4  +  -Si S's |2  (1.A12) 

It  is  possible  to  determine  a  relationship  between  the  remaining  Stokes  con¬ 
stants,  | -Si |,  |<S2|,  and  | -S4 1  from  three  different  boundary  conditions 

l$P  =  j(e"-  -  l)|S4p 

|5.|2  =  yt(e’’-|5<|2-l)  (1.A13) 

In  order  to  make  further  progress,  let  us  consider  the  case 

A  =  1,  B  =  SU  C  =  S2  +  SiS 3,  D  =  0  (1.A14) 

for  which  we  find 

J  =  f  /+-/+  Ixl2  Arg(z)  =  0 

J  1  /_|5,|V^-/_|52  +  5153|V-  Arg(z)  =  n  (LA15) 

where  we  have  defined  X  =  S4  +  S1S5  -f  (S2  +  S1S3).  To  extract  the  value 
of  x,  we  may  consider  the  case  with  D  =  x  with  A,  B  and  C  as  in  (1.A14). 
Because  0+  is  maximally  subdominant  in  the  lower  half  plane  we  may,  in 
fact,  determine  the  value  of  the  solution  at  Arg(z )  =  —n.  The  coefficient  of 
the  maximally  subdominant  solution  must  remain  unchanged  in  the  absence 
of  any  of  the  more  dominant  solutions  so  that 


/+  Arg(z)  =  0 

/+e7?+  Arg(z)  =  —n 

/+e-7?+  +  /—  l-S*!  |2e— T?-  Arg(z)  =  7 r 

-/-l^  +  SrSalV— /+|x|V* 


(1.A16) 


where  we  have  explicitly  kept  the  dependence  of  g+.  Moreover,  we  realize 
that  for  this  situation  A  rather  than  vanishing  takes  a  large  value.  For  the 
case  of  non-zero  A  the  conservation  law  yields 


A 


f+(ev+  -  1)  (1.A17) 

U(l  -  e-r?+  +  |x|V+)  +  /_(|52  +  5i53|V-  -  |-?,|2e-"-) 
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In  the  limit  of  large  t)+ ,  we  realize  that  the  Stokes  constants  which  appear  in 
the  equation  are  finite  so  that 

M2=l  (1.A18) 

Then,  we  determine  from  Eq.  (1.A15)  that 

|S'2  +  51S'3r  =  |51|2e-2’J-  (1.A19) 

Finally,  two  more  boundary  conditions  yield  an  equation  relating  15^1,  IS4I, 
l^l,  IS'sl,  \S2  +  'S’i'5'3|,  and  \s4  +  The  coefficients  |S3|  and  |S’5|  are 

known.  The  coefficients  |52|,  |£4|,  and  |S,2  +  S'iS'3|  may  be  eliminated  in  favor 
of  |Si|.  Then  the  above  equation  along  with  the  physical  Eq.  (1.A10)  is 
sufficient  to  specify  the  coefficients.  We  find  that 

isip  =  is2|2  =  ffy- - 1) 

|54|2  =  1 

\S2  +  SyStf  =  Y~e~r,~(l  —  e~v~) 

l^  +  ^S-sl2  =  e-2r?_  (1.A20) 

The  coefficients  may  be  immediately  determined  from  Eq.  (1.A12) 

Cjc  =  1-e-”-  ,  Rjc  —  e~2r>~  ,  Cfc  =  e~n~(l  —  e~v~)  (1.A21) 


6  Appendix  B:  Equivalency  of  Coupling 
and  Budden  Tunneling  for  Large  n 


d  =  Wo-vW) 

f  =  2(^-1) 

9  =  2((i/2  -  1)k2  +  2t/2.) 


(1.B1) 


(1.B2) 
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and 


^2  ^cl/^c2?Pco  ^cl  /  ^co  ^  —  ^cl/^ii  (l.B3^ 

are  various  mass  and  density  ratios.  We  recognize  this  potential  to  be  a 
Budden  equation  with  a  cutoff  and  resonance.  With  the  appropriate  change 
of  variables 

x  =  “(/  +  9Z)  (1-B4) 

we  recover  the  traditional  form  of  the  Budden  equation  (1.70)  with 

=  d/g3,a  =  cg/d  -  f  (1.B5) 

7  Appendix  C:  Relevant  Frequencies  in  a 
Two-ion  Species  Plasma 

Wave  propagation  in  a  cold,  two-ion  species  plasma  is  characterized  by  the 
gyrofrequencies,  the  cutoff  frequency 


(jjco  —  (niflC2  +  ri2 £}ci)/(fti  ~f  ^2) 

(l.Cl) 

the  crossover  frequency 

Ucr  ~  (nl^c2  +n2^cl)/(nl  +  n2) 

(1.C2) 

and  the  ion-ion  hybrid  frequency 

^01^02(^1  +  w2) 

^ii  ~  Mco 

^1  *  *cl  +  n2*'c2 

(1.C3) 

These  frequencies  have  the  well  defined  ordering 

^cl  ^  <C  Mco  <C  i^cr  ^  ^c2 

(1-C4) 

34 


Part  2 

A  Numerical  Study  of  Mode 
Conversion 

1  Introduction 

In  part  1  of  this  report  we  discussed  the  asymptotic  behavior  of  low  frequency 
modes  which  propagate  at  frequencies  near  the  heavy  ion  gyrofrequency  in  a 
multispecies  plasma.  In  that  work  it  was  our  purpose  to  describe  the  behavior 
of  waves  which  propagate  from  a  region  of  low  magnetic  field  into  a  region 
of  large  magnetic  field.  We  found  that  the  wave  equations  are  amenable  to  a 
Wentzel-Kramers-Brillouin  (WI<B)  analysis  “far”  from  the  crossover-cutoff- 
resonance  triplet  which  occurs  near  the  heavy  ion  gyrofrequency;  however, 
there  exists  a  regime  within  which  such  an  analysis  is  inadequate,  and  the 
propagating  modes  couple  strongly. 

We  characterized  the  behavior  of  various  incident  modes  across  this  region 
by  various  transmission  and  absorption  coefficients  which  corresponded  to  the 
Poynting  flux  attributable  to  each  mode.  From  the  value  of  these  coefficients 
it  is  possible  to  infer  the  power  apportioned  to  each  mode  given  a  wave 
incident  in  a  particular  mode.  In  addition,  we  found  that  substantial  wave 
power  is  absorbed  within  this  region.  One  physical  interpretation  of  such 
absorption  is  that  it  provides  a  source  for  ion  heating  (perhaps  through  mode 
conversion  to  a  thermal  branch  which  can  damp  on  the  ions).  In  addition, 
the  thermal  population  should  also  damp  on  the  transmitted  ion-cyclotron 
wave  providing  an  additional  source  for  heating. 

In  part  1  of  this  report  we  determined  the  primary  transmission  coefficient 
in  all  cases,  and  in  the  case  of  a  large  “gap,”  we  determined  the  reflection 
coefficients  as  well.  In  a  similar,  but  different,  approach  Le  Queau  and  Roux 
[1992]  obtained  the  coefficients  for  the  case  of  a  small  “gap”  and  quantified 
the  associated  absorption.  For  the  case  of  large  angle  of  propagation,  we 
unambiguously  determined  the  coefficients  for  the  ion-cyclotron  mode  (the 
magnetosonic  mode  does  not  propagate  in  that  regime).  In  all  cases,  the 
transmission  is  described  by  Budden-like  transmission  coefficients  and  a  def¬ 
inite  absorption  characterized  the  transmission  process. 
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In  order  to  expand  upon  the  results  which  we  have  obtained  analytically 
and  to  unambiguously  determine  the  wave  power  associated  with  each  mode, 
we  present  a  numerical  solution  to  the  differential  equation  associated  with 
this  problem.  In  particular,  in  the  regime  where  the  tunneling  “gap”  is  small 
we  obtain  the  exact  solutions  to  the  differential  equation  and  from  these 
extract  the  coefficients. 

There  are  a  number  of  advantages  to  be  gained  by  considering  a  numerical 
solution.  Within  the  context  of  a  numerical  solution,  it  is  possible  to  retain 
the  exact  expressions  for  the  Stix  functions  [Sti'x,  1962].  In  order  to  obtain  an 
analytical  solution  we  linearized  these  functions.  Within  that  approximation 
it  is  only  possible  to  retain  one  ion  resonance  and  one  saddle  point.  However, 
as  discussed  in  part  1,  two  saddle  points  exist  at  which  coupling  occurs.  For 
small  or  large  k  one  or  the  other  saddle  point  dominates  the  mode  conversion 
process,  but  for  intermediate  k  both  saddle  points  are  important  and  the 
transition  between  them  should  be  better  understood.  As  we  shall  see,  the 
presence  of  the  extra  terms  does  complicate  the  mode  conversion  process 
slightly.  In  addition,  we  shall  see  that  at  low  altitudes  the  magnetosonic 
and  ion-cyclotron  branches  couple  strongly  (this  is  the  result  of  the  coupling 
point  at  zero  frequency  described  in  the  part  1)  so  that  the  transmitted 
wave  amplitudes  tend  to  be  functions  of  altitude.  Finally,  we  can  model 
the  variations  in  the  magnetic  field  in  a  more  realistic  manner  with  a  cubic 
dependence  on  altitude  rather  than  a  simple  linear  calculation. 

In  the  remainder  of  this  report,  we  develop  and  implement  a  reasonable 
procedure  for  extracting  the  values  of  the  coefficients  numerically.  We  then 
compare  these  numerical  results  with  the  analytical  results  presented  in  part  1 
of  this  report. 


2  Numerical  Procedure 

2.1  Basic  Equations 

In  the  analytic  part  we  considered  a  specific  problem  in  which  waves  propa¬ 
gate  in  the  magnetosphere  from  a  region  of  small  magnetic  field  to  a  region 
of  large  magnetic  field.  As  discussed  in  part  1,  it  is  reasonable  as  a  first 
approximation  to  consider  wave  propagation  in  a  cold  multispecies  plasma. 
Furthermore,  it  is  possible  to  reduce  the  problem  to  one  dimension  if  we 
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consider  variation  only  along  the  direction  of  the  magnetic  field.  Such  an  ap¬ 
proximation  is  reasonable  in  the  specific  context  that  we  considered  because 
parallel  gradients  are  much  larger  than  perpendicular  gradients. 

Given  these  approximations  we  derived  a  fourth-order  coupled  differential 
system  which  appropriately  describes  the  electric  field  components  perpen¬ 
dicular  to  the  magnetic  field 


eV'  +  M#  =  0 


(2.1) 


where 


M  = 


r  —  k2/2  k2/2  \ 

k2/2  1-k2/ 2  ) 


(2.2) 


The  two  components  of  'fr  (tp  and  <f>)  are  the  circularly  polarized  fields  E _ 
and  E+  which  are  functions  of  the  dimensionless  coordinate,  z  =  xj Lb 
where  Lb  is  a  conveniently  chosen  spatial  scale.  We  have  normalized  the 
Stix  functions,  R  and  L,  to  the  large  ratio,  c2 /v2A0  ( vao  is  the  Alfven  velocity 
at  u>  =  flci)>  so  that  the  functions 


r  .1  flci  _ (cj  zb  ujco ) _  f2  3) 

l  9a  uco  (u;  ±  ±  ttc2) 

are  the  order  of  unity.  The  function  gA  contains  the  dependence  of  the  Alfven 
velocity  on  the  coordinate,  z;  f2i,2  are  the  heavy  and  light  ion  gyrofrequencies 
respectively;  and  u >co  =  (n]flC2  +  n2flci)/(ni  +  n2)  is  the  cutoff  frequency. 
The  differential  Eq.  (2.1)  is  characterized  by  two  other  parameters:  the 
small  WI<B  parameter,  t  =  I/IcaoLb,  and  the  normalized  perpendicular 
wavevector,  k  =  k^/lcAo  where  Lao  =  u/vao  is  a  typical  wavevector. 

In  our  approximation,  the  renormalized  Stix  functions  (2.3)  depend  only 
on  the  magnetic  field  which  we  will  define  by  the  function 


g(z)  =  u/ttcl(z) 


(2.4) 


The  particular  form  of  the  Stix  functions 

r  _  v_cojEg_  g 
I  (l  ±  g)(v2  ±  g)  gA 

then  only  depends  on  the  function  g  and  the  obvious  mass  and  density  ratios 
implied  by  the  factors,  u,  which  are  assumed  to  be  constant.  The  function 


(2.5) 
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9a  varies  according  to  the  Alfven  velocity  along  the  field  line  (5n~1/2).  For 
simplicity,  we  shall  assume  throughout  the  analysis  that  g a  remains  constant 
(i.e.  that  the  Alfven  velocity  remains  fixed).  Such  an  assumption  is  certainly 
reasonable  if  the  region  of  coupling  is  localized.  In  addition,  such  an  assump¬ 
tion  also  allows  us  to  match  the  WKB  solutions  more  easily  in  the  region 
of  high  field  (low  altitude).  A  more  relevant  manner  in  which  to  incorpo¬ 
rate  the  actual  decrease  in  the  Alfven  velocity  would  be  to  incorporate  mode 
solutions  on  the  low  altitude  side  of  the  resonance  whose  eigenvalues  are  de¬ 
termined  by  appropriate  ionospheric  boundary  conditions  [ Lysak  and  Dum , 
1983;  Lysak ,  1985].  Although  perhaps  slightly  less  realistic,  our  assumption, 
9a  ^  1?  should  be  sufficient  to  determine  the  coefficients  to  first  order,  and 
it  is  to  be  understood  that  inclusion  of  ionospheric  boundary  conditions  is 
beyond  the  scope  of  the  present  work. 

The  function,  g ,  contains  all  the  essential  information  with  regard  to 
magnetic  field  variations.  The  most  reasonable  dependence  for  g  on  altitude 
is  a  dipole  field.  Thus,  we  expect  that  in  physical  coordinate,  x,  and  scaled 
coordinate,  z,  the  function,  takes  the  form 

LI  1 

9  =  -?  =  7  <2-6> 

where  Lb  is  the  factor  which  appears  in  e.  The  proper  choice  for  Lb  in  this 
instance  is  the  altitude  at  which  the  wave  frequency  matches  the  gyrofre- 
quency  (the  order  of  Re)- 

Having  thus  relegated  all  functional  dependence  into  the  functions  r  and 
l  through  the  function,  g ,  we  solve  the  differential  Eq.  (2.1)  numerically. 
However,  as  discussed  at  length  in  part  1  of  this  report,  the  differential 
equation  has  a  regular  singular  point  at  the  location  of  the  plasma  resonance, 
namely  where  z  =  1.  For  large  values  of  z,  l  >  r,  and  at  the  local  crossover 
frequency  r  =  l.  The  function  l  is  characterized  by  a  zero  at  the  local  cutoff 
frequency  and  a  pole  at  the  heavy  ion  resonance.  The  solution  to  Eq.  (2.1) 
is  multivalued  and  has  a  branch  point  at  z  —  0,  so  that  in  order  to  retain  the 
appropriate  physical  behavior,  we  must  continue  the  solution  above  the  pole 
which  may  be  physically  interpreted  as  damping.  Eq.  (2.1)  must  therefore 
be  suitably  extended  into  the  upper  half  of  the  complex  plane  so  that  the 
pole  may  be  circumvented. 
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2.2  Numerical  Technique 

The  technique  that  we  use  to  solve  Eq.  (2.1)  is  in  principle  straightforward. 
Eq.  (2.1)  constitutes  a  system  of  eight  fourth-order  equations  for  the  real 
and  imaginary  part  of  $ .  One  manner  in  which  to  continue  the  solution,  'k, 
from  the  negative  real  axis  to  the  positive  real  axis  is  to  solve  the  differential 
equation  along  some  parameterized  path  in  the  complex  plane,  z(r),  char¬ 
acterized  by  a  real  parameter  r.  The  solution  is  then  obtained  by  solving 
the  equivalent  eight  first-order  ODE’s  in  the  parameter  r  by  means  of  a  well 
known  integrating  scheme  such  as  the  Runge-Kutta,  Bulirsch-Stoer,  or  the 
predictor-corrector  method  [Press  et  al .,  1986].  The  generalization  of  Eq. 
(2.1)  to  a  continuous  (although  not  necessarily  smooth)  path  is  straightfor¬ 
ward  and  is  discussed  in  detail  elsewhere  [Johnson,  1992]. 

In  order  to  obtain  the  correct  branch  of  the  multivalued  solution  of  Eq. 
(2.1)  we  must  continue  the  solution  above  the  pole  at  z  =  0  in  the  complex 
plane.  Therefore,  we  solve  the  differential  equation  along  a  parameterized 
path  which  passes  above  the  pole.  Suppose  that  we  know  the  solution,  ty,  at 
some  point,  zq,  along  the  negative  real  axis.  In  order  to  determine  the  value 
of  the  solution  at  some  point,  z ,  along  the  real  axis,  we  may  integrate  the 
equivalent  differential  equation  along  the  parameterized  path  such  as  that 
illustrated  in  Figure  7. 

z\(t)  =  t  z0  <  t  <  1  —  p 

z2(t)  -  pex  p  ^(1  +  1  -  T)^  l-p<T<l+p  (2.7) 

z3 (t)  =  t  1  +  p  <  r  <  z 

which  consists  of  integration  along  the  real  axis  from  z0  to  1  —  p  (recall  that 
the  pole  is  at  z  =  1)  followed  by  integration  along  a  circular  arc  above  the 
pole  to  1  +p  on  the  real  axis.  Finally,  the  solution  at  z  may  be  obtained  from 
further  integration  along  the  real  axis.  The  solution  must  not  depend  on  the 
choice  of  p  or  for  that  matter  upon  the  choice  of  parameterized  path.  This 
requirement  provides  an  additional  check  on  the  reliability  of  the  procedure. 

A  general  solution  to  the  differential  equation  may  be  obtained  by  taking 
four  different  initial  conditions  at  the  point  Zo  which  correspond  to  four  lin¬ 
early  independent  solutions.  Once  we  have  found  a  general  solution  to  the 
differential  equation  we  may  decompose  it  into  its  constituent  WKB  solu¬ 
tions.  Far  from  the  coupling  region,  the  WKB  solutions  are  asymptotically 
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Figure  7:  The  parameterized  path  along  which  we  solve  the  differential 
equation  consists  of  three  contours.  The  first  integrates  along  the  real  axis 
to  some  point  p.  To  continue  the  solution  to  the  positive  real  axis,  we  take 
a  contour  along  a  circular  arc  above  the  pole.  Finally  the  solution  may  be 
obtained  by  continuing  the  solution  along  the  real  axis. 


close  (in  the  parameter  e)  to  the  actual  solutions. 


2.3  Extracting  the  Coefficients 

Once  four  linearly  independent  solutions  for  the  lower  component,  0,  of  $ 
have  been  calculated  (xpj,  j  =  1,4)  using  the  method  detailed  in  the  last 
section,  it  is  relatively  straightforward  to  determine  the  coefficients.  In  the 
asymptotic  regime,  the  WKB  analysis  yields  four  solutions  (0j,  j  =  1,4). 
In  the  asymptotic  region  far  to  the  left  of  the  coupling  region  along  the 
negative  real  axis  and  far  to  the  right  along  the  positive  real  axis  the  WKB 
solutions  are  reasonable  approximations  to  the  actual  solutions  and  serve  as 
basis  functions  as  illustrated  in  Figure  8. 

The  general  solution  to  the  upper  component  of  (2.1)  may  be  written  in 
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Figure  8:  In  regions  far  from  the  pole  the  WKB  solutions  are  asymptoti¬ 
cally  valid.  The  numerical  solution  connects  the  regions  in  which  the  WKB 
solutions  are  accurate. 


the  form 

Tp  =  CXjlpj 


(2.8) 


where  we  have  introduced  the  summation  notation.  Four  appropriate  bound¬ 
ary  conditions  are  sufficient  to  determine  a  particular  solution  of  the  differen¬ 
tial  equation.  The  boundary  conditions  which  we  impose  on  the  solution  are 
related  to  the  coefficients  of  the  WKB  solutions  above  and  below  the  coupling 
region.  Because  the  WKB  solutions  are  basis  functions  in  the  asymptotic 
regime,  the  solution,  ip,  may  be  decomposed  as 


ip  =  7* 0j  ,  -»  ±oo 


(2.9) 


where  the  decomposition  may  be  different  along  the  positive  and  negative 
real  axis  in  light  of  the  Stokes  phenomenon  [ Heading ,  1962].  In  the  following 
analysis  we  will  refer  to  the  region  3?(^)  — >  +oo  as  “above”  (low  field  side) 
and  we  will  describe  waves  propagating  from  3?(^)  >  0  to  3?(z)  <  0  as 
“downgoing”  waves  (“below”  and  “upgoing”  have  the  obvious  meanings). 
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The  boundary  conditions  are  related  to  the  power  flow  through  the  con¬ 
served  quantity 

J  =  t(<Pt'tf  -  (2.10) 

i.e.,  the  field  aligned  Poynting  flux,  described  in  part  1.  As  detailed  in  Ap¬ 
pendix  A,  the  coefficients,  7,  correspond  to  power  flux  in  the  four  propagating 
modes.  To  clearly  specify  the  boundary  conditions  we  choose  Qj  as  speci¬ 
fied  in  Appendix  A  corresponding  for  parallel  propagation  to  a  downgoing 
left-hand  circularly  polarized  (LHCP)  wave,  upgoing  LHCP  wave,  downgoing 
right-hand  circularly  polarized  (RHCP)  wave,  and  an  upgoing  RHCP  wave 
respectively. 

The  values  of  7^  must  be  consistent  with  the  boundary  conditions  related 
to  our  problem.  For  a  downgoing  magnetosonic  wave  we  require  that 


0  ) 

c,  ^ 

C2 

0 

1 

5 

c3 

C4  / 

f 

l  0  J 

(2.11) 


Physically,  these  conditions  correspond  to  no  upgoing  waves  from  below  the 
coupling  region,  and  above  the  coupling  region,  the  only  incident  wave  is 
the  downcoming  approximately  RHCP  mode.  The  amplitudes,  cjt  are  to  be 
determined  and  are  related  to  the  power  transmitted  in  the  incident  mode, 
Tms  (c3);  reflected  in  the  incident  mode,  RMS  (c4);  transmitted  in  the  coupled 
mode,  Cj^s  (c2);  and  reflected  in  the  coupled  mode,  C^5  (ci). 

For  an  incident  ion-cyclotron  wave,  we  require  that  the  amplitudes,  7±, 
take  the  form 


1  ^ 

4  \ 

4 

0 

0 

C3 

<  ) 

'v  0  J 

consonant  with  the  physical  boundary  conditions.  In  this  case  c'j  correspond 
to  transmitted,  Tjc ;  reflected,  Rjc ;  coupled  transmitted,  Cjc ;  and  coupled 
reflected,  C'pc ,  wave  power  respectively. 

In  order  to  determine  the  coefficients,  aj,  we  decompose  the  linearly  in¬ 
dependent  functions,  ipj,  into  the  constituent  WKB  solutions  on  both  the 
positive  and  negative  imaginary  axis.  In  light  of  the  Stokes  phenomena,  the 
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expansions  are  not  the  same.  Then 

rpj  =  fi%Qk  ,  »(*)^±<x>  (2.13) 

which  determines  the  relationship  between  the  coefficients 

it  =  (2-14) 

The  Eq.  (2.14)  constitute  a  system  of  8  complex  equations  in  8  unknown 
complex  variables.  The  variables  which  must  be  determined  are  the  values  of 
o.j  and  the  amplitudes,  cj  or  c' .  The  values  of  a  are  obtained  from  the  four 
equations  with  definite  boundary  conditions.  The  values  of  the  amplitudes 
are  determined  by  the  remaining  four  equations.  The  amplitudes  determine 
the  power  flow  and,  thus,  the  coefficients  T,  R  CT ,  and  CR  as  detailed  in 
Appendix  A.  In  either  the  case  of  incident  magnetosonic  or  ion  cyclotron 
waves,  the  absorption  coefficient  is  of  the  form 

A  =  1  ~(T+R  +  Cr  +  Ct)  (2.15) 

The  physical  interpretation  of  the  coefficients  is  discussed  in  Figures  1.3  and 
1.4. 

3  Numerical  Complications 

In  studying  the  ideal  system  of  Eq.  (2.1)  with  the  well  approximated  spatial 
dependence  on  the  magnetic  field,  it  is  important  for  us  to  address  some  basic 
limitations  in  the  numerical  analysis.  Numerically,  the  equations  are  suffi¬ 
ciently  tractable  that  the  exact  solution  to  the  well  approximated  Eq.  (2.1) 
may  be  determined  to  arbitrary  accuracy  (within  the  limits  of  computer  ac¬ 
curacy).  We  can  make  a  statement  such  as  this  for  a  number  of  reasons. 
We  have  solved  the  differential  equation  using  three  different  techniques, 
Bulirsh-Stoer,  Runge-Kutta  and  predictor-corrector  [Press  et  a/.,  1986].  All 
techniques  produce  the  same  results  to  within  the  limitations  of  double  preci¬ 
sion  accuracy.  Moreover,  we  have  transformed  the  Eq.  (2.1),  using  a  Ricatti 
Transform,  into  a  set  of  first-order  non-linear  differential  equations  [Smith 
and  Whitson ,  1978].  (Such  a  transformation  is  particularly  useful  for  sup¬ 
pressing  exponentially  growing  solutions.)  The  numerical  solution  that  we 
have  obtained  to  the  transformed  Ricatti  equation  is  again  consistent  with 
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the  numerical  solution  of  (2.1)  to  arbitrary  accuracy.  In  principle,  then,  the 
accuracy  of  the  coefficients  that  we  determine  is  limited  to  the  accuracy  of 
the  imposed  boundary  conditions. 

The  boundary  conditions  which  we  impose  upon  the  numerical  solution 
are  at  best  limited  to  the  validity  of  the  WKB  solutions.  That  is,  the  accu¬ 
racy  of  the  coefficients  is  limited  by  the  accuracy  of  p±k.  Once  a  numerical 
solution  xfjj  is  obtained,  the  most  straightforward  manner  in  which  to  deter¬ 
mine  the  coefficients,  is  to  compare  the  value  of  the  numerical  solution, 

with  the  four  WKB  solutions,  0^,  at  four  different  locations.  Alterna¬ 
tively,  we  could  match  i\)j  and  its  derivatives  at  a  single  point  with  the  0^  and 
their  derivatives.  More  elaborate  schemes  are  also  conceivable  in  which  xj)j  is 
matched  in  a  number  of  locations  and  an  average  taken.  Moreover,  because 
we  consider  only  the  leading  order  WKB  approximation,  it  is  important  that 
we  match  the  numerical  and  WKB  solutions  in  a  region  small  compared  to 
a  typical  wavelength.  On  the  other  hand,  if  this  region  is  taken  too  small, 
then  the  values  of  0^  do  not  vary  sufficiently  to  determine  the  coefficients 
Vjh  to  the  required  accuracy.  Indeed,  the  quality  of  the  matching  is  a  basic 
limitation  of  our  analysis.  In  the  results  which  follow,  we  have  matched  the 
solutions  in  a  region  such  that  the  WKB  approximation  is  reasonably  valid, 
and  we  have  varied  the  width  of  the  matching  region  and  found  that  within 
the  upper  and  lower  limits  mentioned  above  deviations  from  the  presented 
results  are  not  substantial. 

For  large  values  of  k,  the  RHCP  magnetosonic  mode  has  very  small  so 
that  matching  to  the  WKB  solutions  becomes  somewhat  suspect.  One  should 
therefore  question  the  exact  manner  in  which  the  coefficients  behave  near  the 
critical  k;,  which  corresponds  to  the  transition  from  four  propagating  WKB 
modes  to  two  propagating  modes,  because  the  meaning  of  matching  to  WKB 
boundary  conditions  is  questionable.  However,  we  should  point  out  that 
the  coefficients  do  behave  reasonably  smoothly  and  it  is  not  unreasonable  to 
assume  that  the  asymptotic  solutions  are  reasonably  accurate  even  when  the 
WKB  smallness  parameter,  e  ~  1  [Bender  and  Orszag ,  1978], 

We  should  also  mention  the  dependence  of  the  solutions  on  the  parameter, 
p,  which  was  introduced  to  circumvent  the  pole  at  z  =  1.  The  solutions 
i’j  have  virtually  no  dependence  on  the  parameter  p  as  long  as  p  is  taken 
sufficiently  small.  In  the  case  of  very  large  p,  the  exponentially  growing 
solution  in  the  upper  half  plane  dominates  the  subdominant  solution  to  the 
machine  precision  in  which  case  an  inaccurate  numerical  solution  is  obtained. 
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For  all  results  presented,  the  value  of  p  is  taken  sufficiently  small  so  that 
numerical  stability  is  retained. 

Finally,  as  2  ->  0  so  that  g  -4  00  we  find  that  r  =  l  which  corresponds 
to  the  lower  crossover  frequency.  Coupling  associated  with  this  crossover 
frequency  occurs  between  2  =  0  and  2  =  1.  Moreover,  the  strong  depen¬ 
dence  of  the  magnetic  field  on  z~3  means  that  this  coupling  occurs  over  an 
extended  range  of  2  (r  «  /  over  a  large  range).  As  a  result,  the  numerical 
values  of  the  coupling  coefficients  are  functions  of  altitude.  However,  the 
physical  boundary  conditions  which  we  impose  on  the  upgoing  waves  below 
the  heavy  ion  resonance  are  independent  of  the  location  at  which  we  match 
the  solutions.  Then  it  turns  out  that  the  reflection  and  absorption  coeffi¬ 
cients  do  not  depend  upon  the  location  at  which  the  solution  is  matched  for 
2  <  1;  however,  the  relative  value  of  the  transmission  coefficients  will  depend 
upon  the  location  at  which  they  are  evaluated  to  a  slight  degree  because  the 
two  downgoing  Alfven  modes  couple  at  low  altitudes.  In  the  following  anal¬ 
ysis  we  match  to  the  slowly  varying  WKB  solutions  at  some  point  between 
2  =  0  and  2  =  1.  In  all  cases,  we  have  chosen  to  evaluate  the  coefficients 
in  a  range  over  which  the  variation  of  the  coupling  coefficients  with  altitude 
is  small.  Indeed,  this  range  is  identically  the  altitude  range  over  which  the 
boundary  conditions  may  be  optimally  described  in  terms  of  the  coefficients 
of  the  WKB  solutions. 

4  Results 

4.1  Parameter  Regime 

In  the  following  sections  we  determine  the  coefficients  for  several  cases  of  inci¬ 
dence  as  described  previously.  The  parameters  which  are  of  greatest  interest 
are  the  magnetic  field,  magnetic  field  scale  length,  and  the  masses  and  den¬ 
sities  of  the  various  constituent  ions.  The  masses  are  fixed  parameters  which 
we  take  to  be  hydrogen  and  oxygen  (both  singly  ionized).  The  magnetic  field 
and  densities  are  embedded  in  the  WKB  parameter  e  and  the  functions  r 
and  l  depend  only  upon  the  density  ratios.  Hence,  a  particular  solution  is 
determined  from  the  specification  of  e  and  the  concentration  of  the  minority 
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species  r).  Typical  values  for  the  parameter  c  are 


e  =  Tl-  ~° .5r-lp~l  (2.16) 

kALB 

where  r  is  the  scale  length  in  earth  radii  (~  1  —  3)  and  p  is  the  mass  density 
of  ions  (~  1  ~  100)  [Moore  et  al. ,  1987,  and  references  therein].  We  will 
explore  the  parameter  space  for  t  ranging  from  approximately  0.2  to  0.03.  In 
this  range,  the  WKB  approximation  is  good  and  should  provide  reasonable 
matching  conditions  to  the  WKB  solutions.  We  will  in  all  cases  concentrate 
on  that  range  of  oxygen  concentration  consistent  with  reasonable  absorption 
levels.  We  have  found  this  range  to  vary  from  0.5  to  5  percent.  In  order  to 
consider  the  limiting  case  of  a  large  “gap”  we  will  take  e  to  be  small  and  rj 
to  be  large. 

4.2  Incident  Magnetosonic  Waves 

For  the  results  of  this  section,  we  impose  the  boundary  conditions  for  down¬ 
coming  magnetosonic  waves  and  obtain  the  coefficients  (2.A7).  At  this  junc¬ 
ture  it  would  be  useful  for  the  reader  to  examine  Figure  3  to  understand 
the  physical  meaning  of  the  coefficients.  In  that  figure  a  wave  is  incident 
on  the  magnetosonic  branch  with  unit  flux.  The  fractions  of  wave  flux  Tms 
and  Rms  are  transmitted  and  reflected  on  the  magnetosonic  branch  while  the 
fractions  C^s  an<^  C ms  are  transmitted  and  reflected  along  the  ion-cyclotron 
branches  as  indicated. 

We  have  obtained  full  analytical  solutions  for  the  limit  of  a  large  “gap.” 
In  Figure  9,  we  show  the  numerical  solutions  to  the  differential  equation. 
Superimposed  on  the  graph  we  have  shown  the  analytical  results  (1.44)  which 
we  have  obtained  from  a  phase  integral  analysis  consistent  with  a  large  “gap” 
for  which  there  is  no  absorption.  The  results  are  in  very  good  agreement.  The 
most  striking  difference  is  that  the  coefficient  Cj,jS  is  non-zero  in  contrast  to 
the  analytic  situation.  However,  the  analytical  results  only  incorporate  the 
existence  of  the  coupling  point  near  the  crossover  frequency.  As  discussed 
in  part  1,  coupling  between  the  two  Alfven  waves  can  be  very  strong  at 
frequencies  below  the  heavy  ion-cyclotron  frequency.  It  is  this  coupling  which 
gives  rise  to  the  wave  power  found  on  the  ion  cyclotron  branch. 

In  Figure  10,  we  compare  the  numerical  transmission  coefficient  for  a 
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Figure  9:  For  the  limit  of  large  “gap”  the  numerical  solutions  and  the 
analytical  results  are  in  reasonable  agreement.  We  have  taken  e=0.03  and 
7 i  =  5%  for  which  there  is  no  absorption.  The  numerical  solutions  are  plotted 
in  heavy  set  type  and  the  analytical  results  are  indicated  in  light  type.  Note 
that  for  the  numerical  results  the  coefficient  exp (— 7r/?+/e)  although  small  is 
nonzero  in  contrast  to  the  analytical  results.  This  behavior  is  the  result  of 
coupling  between  the  two  Alfven  waves  below  the  ion  resonance  correspond¬ 
ing  to  the  coupling  point  at  zero  frequency. 

small  “gap”  with  the  analytical  result  (1.44).  We  have  plotted  the  numer¬ 
ical  transmission  coefficient  in  heavy  solid  and  the  analytical  coefficient  in 
light  solid.  The  analytical  results  predict  no  coupled  wave  Cj^s.  However, 
coupling  between  the  two  modes  does  occur  below  the  heavy  ion  resonance. 
Because  this  coupling  does  not  involve  reflection  one  might  expect  the  total 
transmission  Tms  +  CmS  to  match  the  analytical  transmission  coefficient.  We 
have  plotted  this  sum  as  a  dashed  line.  Indeed  the  sum  compares  well  with 
the  analytical  results. 

In  the  regime  which  corresponds  to  a  large  “gap”  our  numerical  results  are 
in  good  agreement  with  the  analytical  results.  Moreover,  even  in  the  limit  of 
a  small  “gap,”  the  analytical  results  are  in  reasonable  agreement.  Let  us  now 
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Figure  10:  Transmission  coefficient  for  the  magnetosonic  wave  in  the  limit 
of  small  “gap”  with  e=0.2  and  77  =  1%.  Note  that  the  total  transmission  is 
in  good  agreement  with  analytical  result. 

explore  the  parameter  space.  In  Figure  1 1  we  have  plotted  the  coefficients 
for  the  case  t  =  0.2  and  77  =  1%.  For  parallel  propagation,  no  coupling 
occurs.  As  the  k  increases  strong  coupling  between  the  two  transmitted 
Alfven  waves  occurs,  and  the  sum  of  of  those  coefficients  reproduces  the 
analytical  transmission  coefficient.  For  large  k  waves  are  reflected  in  the  two 
upward  propagating  modes.  The  cutoff  frequency  increases  substantially  for 
large  enough  k  and  the  “gap”  becomes  too  large  to  penetrate.  In  this  case  the 
reflection  coefficient  grows  rapidly,  and  for  k  ~  0.95  the  magnetosonic  wave 
no  longer  propagates  (the  wave  is  reflected  at  the  lower  hybrid  frequency 
well  above  the  hydrogen  gyrofrequency).  For  a  substantial  range  of  k  strong 
absorption  occurs  peaking  at  20%  for  large  k. 

For  smaller  densities,  the  absorption  increases  to  about  25%  although  it 
tends  to  sharpen  and  occur  mostly  at  large  k.  Indeed,  as  the  minority  species 
density  tends  to  zero,  the  transmission  coefficient,  Tms ,  is  approximately  1 
and  all  other  coefficients  vanish.  This  means,  in  effect,  that  for  very  small 
oxygen  densities  the  waves  propagate  as  if  there  were  no  oxygen  present. 
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Figure  1 1 :  Coefficients  for  the  case  e  =  0.2  and  rj  =  1%. 


For  larger  densities  (77  =  5%),  as  illustrated  in  Figure  12,  we  find  that  the 
absorption  diminishes  rapidly  although  the  waves  still  couple  at  low  altitude. 
In  addition,  we  find  that  total  reflection  occurs  at  a  much  smaller  value  for 
k.  Absorption  appears  to  be  very  small.  However,  one  should  keep  in  mind 
that  in  the  case  of  a  large  gap  all  power  upcoming  in  the  LHCP  mode  will 
be  absorbed.  If  perfect  reflection  were  to  occur  at  ionospheric  altitudes,  it  is 
conceivable  that  a  substantial  portion  of  C^fS  would  also  be  absorbed. 

For  smaller  values  of  e  (=0.08),  as  illustrated  in  Figure  13,  we  find  that 
less  coupling  occurs  between  the  two  downgoing  transmitted  waves.  In  addi¬ 
tion,  reflection  is  much  stronger  at  intermediate  values  of  k.  This  is  because 
the  “gap”  contains  more  wavelengths  and  transmission  is  substantially  di¬ 
minished.  Reasonably  strong  absorption  still  occurs  over  a  vast  range  of  k 
and  is  typically  about  7  percent.  It  should  also  be  noted  that  a  substantial 
portion  of  the  reflected  wave  Cj^s  can  also  contribute  to  heating  upon  re¬ 
flection.  For  even  smaller  values  of  the  parameter  t  (=0.03)  we  find  that  the 
absorption  becomes  negligible  for  rj  >  1%,  however,  the  absorption  remains 
reasonably  significant  for  77  =  0.5%. 

In  summary,  we  have  found  that  incident  magnetosonic  waves  couple 
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Figure  12:  Coefficients  for  the  case  t  =  0.2  and  rj  =  5%. 


strongly  near  the  crossover  frequency  for  reasonable  values  of  the  minority 
species  density.  For  larger  values  of  ft  the  waves  are  completely  reflected,  and 
strong  absorption  occurs  over  a  substantial  range  of  ft. 

4.3  Incident  Ion- Cyclotron  Waves 

For  the  results  of  this  section,  we  impose  the  boundary  conditions  for  down¬ 
coming  magnetosonic  waves  and  obtain  the  coefficients  (2.A8).  For  values 
of  k  smaller  than  the  critical  ft  for  which  the  magnetosonic  mode  is  cut  off 
we  determine  these  values  numerically.  For  larger  values  of  ft  we  use  the 
analytical  results  obtained  in  part  1.  Numerical  extraction  of  the  coefficients 
in  the  case  of  large  ft  involves  substantial  complications.  In  this  section  we 
will  obtain  correct  estimates  of  the  absorption  for  both  small  and  large  val¬ 
ues  of  ft,  and  we  will  discuss  the  meaning  and  correctness  of  evaluating  the 
coefficients  for  intermediate  ft. 

First,  let  us  consider  the  coefficients  for  values  of  ft  smaller  than  the  criti¬ 
cal  ft  for  which  the  magnetosonic  mode  is  cut-off.  In  order  to  understand  the 
meaning  of  these  coefficients  is  illustrated  in  Figure  4.  A  wave  is  incident  on 
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Figure  13:  Coefficients  for  the  case  c=0.08  and  r]  =  1%. 

the  ion-cyclotron  branch  with  unit  flux.  The  fractions  of  wave  flux  Tic  and 
Ric  are  transmitted  and  reflected  on  the  ion-cyclotron  branches  while  the 
fractions  Cfc  and  Cjc  are  transmitted  and  reflected  along  the  magnetosonic 
branch  as  indicated. 

We  have  obtained  full  analytical  solutions  for  the  limit  of  a  large  “gap.” 
In  Figure  14  we  show  the  numerical  solutions  to  the  differential  equation. 
Superimposed  on  the  graph  we  have  shown  the  analytical  results  (1.43)  which 
we  have  obtained  from  a  phase  integral  analysis  consistent  with  a  large  “gap” 
for  which  there  is  no  absorption.  As  for  the  case  of  an  incident  magnetosonic 
wave,  the  results  are  in  very  good  agreement.  Again  a  small  amount  of  wave 
flux  is  found  in  the  coupled  transmitted  wave. 

In  Figure  15  we  compare  the  analytic  solution  for  a  small  “gap”  with 
the  numerical  calculation.  Once  again  coupling  between  the  two  transmitted 
modes  diminishes  the  power  found  on  the  transmitted  ion-cyclotron  mode. 
However,  the  results  are  certainly  quantitatively  correct  for  small  values  of 
k  in  which  case  the  low  frequency  coupling  point  does  not  enter  into  the 
calculation.  For  larger  the  results  are  qualitatively  correct  and  are  good 
for  an  order  of  magnitude  estimate;  however,  the  existence  of  two  different 
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Figure  14:  For  the  limit  of  large  “gap”  the  numerical  solutions  and  the 
analytical  results  are  in  reasonable  agreement.  We  have  taken  e  =  0.03  and 
r]  =  5%  for  which  there  is  no  absorption.  The  numerical  solutions  are  plotted 
in  heavy  set  type  and  the  analytical  results  are  indicated  in  light  type. 


coupling  points  (low  frequency  and  crossover — see  Figure  2)  cannot  be  in¬ 
corporated  into  a  linear  model  and  the  analytic  results  are  most  accurate 
when  the  coupling  is  dominated  by  one  of  the  two  coupling  points. 

In  Figure  16  we  have  plotted  the  coefficients  for  the  case  e  =  0.2  and 
t?  =  1%.  For  parallel  propagation,  the  problem  is  described  as  a  Budden 
tunneling  problem  for  the  ion-cyclotron  branch,  and  the  ion-cyclotron  and 
magnetosonic  branches  do  not  couple.  For  larger  angles  of  propagation,  the 
tunneling  increases  as,  in  fact,  predicted  by  the  analytical  analysis  in  Fig¬ 
ure  15.  For  larger  values  of  k,  coupling  between  the  downgoing  modes  occurs 
at  the  crossover  frequency,  and  substantial  wave  energy  is  transferred  to  the 
coupled  magnetosonic  wave.  Indeed,  for  the  case  of  a  small  “gap”  as  in 
Figure  14,  coupling  to  the  coupled  transmitted  wave,  CJci  is  complete  for 
large  values  of  k.  Substantial  absorption  occurs  over  a  very  large  range  of  k 
taking  a  value  over  20%  over  most  of  the  range.  In  Figure  17  we  explore 
the  transmission  coefficients  for  parallel  propagation  for  c  =  0.2  (for  parallel 
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Figure  ,15:  Transmission  coefficient  for  the  ion-cyclotron  wave  in  the  limit  of 
small  “gap”  with  e  =  0.2  and  77  =  5%.  Note  that  the  transmission  is  in  best 
agreement  for  small  k. 

propagation  the  two  modes  are  uncoupled  so  that  Cfc  =  Cjc  —  0)-  As  is 
clear,  substantial  absorption  occurs  over  a  large  range  of  densities  ranging 
from  77  =  0.5%  to  77  =  4%.  For  larger  values  of  k,  the  absorption  remains 
essentially  the  same  as  for  parallel  propagation  up  to  reasonably  large  values 
of  K. 

As  for  the  magnetosonic  case,  substantial  absorption  still  occurs  for  e  = 
0.08  as  indicated  in  Figure  18.  In  this  case,  the  “gap”  is  large  so  that 
most  of  the  wave  power  is  reflected  at  small  angles  of  propagation  rather 
than  transmitted;  however,  for  larger  angles  of  propagation,  tunneling  in 
both  modes  increases  substantially  until  most  of  the  power  is  coupled  into 
the  magnetosonic  branch.  For  the  case  of  parallel  propagation  with  e  =  0.08, 
the  peak  absorption  levels  are  skewed  toward  smaller  values  of  density  (77  = 
0.5%).  For  smaller  values  of  e  the  absorption  is  completely  negligible  and  the 
wave  is  completely  reflected. 

For  large  values  of  k,  we  have  shown  using  a  variety  of  methods  that 
the  transmission  properties  are  described  as  a  Budden  tunneling  problem 
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Figure  16:  Coefficients  for  the  case  t  =  0.2  and  rj  ~  1%. 


at  the  ion-ion  hybrid  frequency.  Let  us  consider  the  analytical  results.  In 
Figure  19  we  have  plotted  in  bold  type  the  transmission  factors  found  in 
Eq.  (1.54,  1.58)  along  with  the  corresponding  absorption  factor.  Physically, 
the  coefficients  correspond  to  Figure  6.  In  light  type  we  have  plotted  the 
transmission  factors  obtained  from  the  embedded  Eq.  (1.71).  They  are  in 
good  agreement  in  the  regime  in  which  they  are  both  valid.  It  is  to  be 
noted  that  there  is  reasonable  agreement  between  the  transmission  factor, 
T/c,  near  the  upper  bound  of  the  numerical  solution  of  Figure  16  and  the 
large  k  solutions  of  Figure  19.  The  absorption  and  reflection  coefficients  do 
not  match  well.  While  it  is  to  be  expected  that  the  dominant  WKB  solution 
would  be  continuous  in  k  it  is  not  necessarily  expected  that  the  coefficient  of 
the  subdominant  solutions  would  remain  continuous. 

For  large  k  the  transmission  coefficients  approach  a  constant  value  deter¬ 
mined  from  Budden  tunneling  at  the  ion-ion  hybrid  frequency.  In  Figure  20 
we  plot  the  coefficients  against  the  concentration  of  oxygen.  As  for  the  case 
of  parallel  propagation  strong  absorption  occurs  for  a  large  range  of  densi¬ 
ties  ranging  from  rj  =  0.5%  to  rj  =  7%.  In  all  cases  the  absorption  remains 
essentially  constant  for  k  >  1. 
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Figure  17:  Coefficients  for  parallel  propagation  with  e  =  0.2.  Strong  ab¬ 
sorption  occurs  for  0.5%  <  rj  <  4%. 


For  smaller  values  of  c  the  absorption  still  remains  very  large  for  rea¬ 
sonable  densities.  As  illustrated  in  Figure  21  although  the  reflection  is 
now  much  larger  than  the  transmission,  the  absorption  remains  over  20fac- 
tor,  Tic ,  again  matches  reasonably  well  with  that  coefficient  for  small  k  (see 
Figure  18),  but  the  other  coefficients  are  again  discontinuous.  From  Fig¬ 
ure  22,  it  is  clear  why  the  transmission  has  not  diminished  from  the  case 
with  e  =  0.2.  The  absorption  is  strongly  peaked  around  t]  =  1.  For  smaller 
values  of  e  the  region  of  strong  absorption  moves  to  very  small  (probably 
unphysical)  values  of  r /. 

In  summary,  we  have  found  that  incident  ion-cyclotron  waves  undergo 
Budden-like  tunneling  at  the  cutoff  frequency  for  nearly  parallel  propagation 
and  at  the  ion-ion  hybrid  frequency  for  nearly  perpendicular  propagation. 
The  absorption  is  very  strong  over  a  large  range  of  n  and  for  a  reason¬ 
able  range  of  eta.  For  intermediate  values  of  k,  coupling  near  the  crossover 
frequency  becomes  important  and  a  substantial  amount  of  wave  power  is 
transmitted  in  the  magnetosonic  mode.  However,  for  large  enough  /c,  the 
magnetosonic  mode  is  cut-off  and  that  wave  power  is  reflected  back  into  the 
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Figure  18:  Coefficients  for  the  case  t  =  0.08  and  77  =  1%. 


system. 

4.4  Validity  of  Analytic  Approximation 

The  analytic  results  that  were  obtained  in  part  1  and  those  obtained  by  Le 
Queau  and  Roux  [1992]  are  in  excellent  agreement  for  small  and  large  k  for 
which  one  or  the  other  of  the  two  saddle  points  dominate  the  mode  conversion 
process.  However,  for  a  small  “gap”  there  is  substantial  coupling  between  the 
two  downgoing  modes  resulting  from  the  cutoff  saddle  point.  The  numerical 
results  have  the  advantage  of  being  able  to  retain  both  saddle  points.  It 
is  also  important  to  point  out  that  the  absorptions  published  by  Le  Queau 
and  Roux  [1992]  are  somewhat  misleading.  They  show  strong  absorption  for 
77  ~  10  —50%.  First,  their  approximation  is  only  valid  for  densities  less  than 
3%.  Secondly,  our  results  clearly  show  that  absorption  is  at  best  limited 
to  densities  below  about  5%.  The  reason  behind  this  discrepancy  is  that 
they  have  fixed  the  ratio  their  tunneling  factor,  rjLB/ A,  so  that  the  WI<B 
approximation  is  actually  not  valid  for  r/  more  than  about  1%. 
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Figure  19:  Analytical  approximations  for  the  transmission  coefficients  for 
t  =  0.2  and  77  =  1%  at  large  k  (which  corresponds  to  large  angles  of  propa¬ 
gation).  We  have  plotted  the  coefficients  related  to  the  phase  integral  trans¬ 
mission,  exp(—7r /?+/e),  in  dark  type.  The  transmission  coefficients  found 
from  the  embedded  Budden  equation  for  tunneling  through  the  “gap”  are 
indicated  in  light  type. 

4.5  Physical  Considerations 

We  have  found  that  the  analysis  of  an  incident  ion-cyclotron  wave  is  some¬ 
what  more  complicated  than  the  magnetosonic  case  in  that  it  involves  an 
interesting  change  in  the  behavior  for  the  coefficients  for  small  and  large  k. 
In  particular,  it  is  important  to  realize  that  the  continuity  of  the  the  trans¬ 
mission  properties  for  these  waves  is  somewhat  dependent  upon  the  model 
from  which  boundary  conditions  are  obtained.  The  fact  that  two  different 
regimes  exist  for  k  <  1  and  k,  >  1  is  represented  by  the  transition  of  the 
cutoff  frequency  where  a  cutoff  appears  below  the  ion-cyclotron  frequency 
(instead  of  above).  Then  the  boundary  conditions  at  ±00  are  determined  by 
two  WKB  solutions  rather  than  four.  One  expects  that  such  a  sudden  tran¬ 
sition  in  the  behavior  of  the  modes  is  not  physically  meaningful  although  it 
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Figure  20:  Coefficients  for  large  k,  with  c  =  0.2.  Strong  absorption  occurs 
for  0.5%  <  r]  <  5%. 

is  certainly  mathematically  feasible  (recall  the  discussion  of  the  sudden  jump 
in  the  Stokes  coefficients  across  a  definite  value  for  Arg(z)). 

In  light  of  our  present  model,  simply  from  physical  considerations  we 
would  expect  to  find  a  discontinuous  jump  in  the  absorption  coefficient  as  the 
cutoff  moves  from  +oo  to  —  oo.  If  the  magnetosonic  mode  propagates  at  —  oo 
then  one  expects  to  find  that  all  the  wave  power  on  the  magnetosonic  mode 
is  transmitted;  however,  if  that  mode  is  reflected  above  the  region  where  the 
coefficients  are  evaluated,  then  all  of  the  wave  power,  Cfc ,  on  that  branch  is 
reflected  back  into  the  system.  Because  Cfc  increases  substantially  at  large 
k  where  coupling  at  the  crossover  frequency  is  most  complete,  this  amount  of 
reflected  energy  is  substantial.  Such  reflection  greatly  modifies  the  reflection 
and  absorption  coefficients  although  it  should  not  affect  the  transmission 
coefficient.  As  we  have  seen,  the  transmission  coefficient  actually  matches 
reasonably  well  across  the  critical  k. 

Physically,  we  might  expect  more  continuous  behavior  for  the  coefficients; 
however,  producing  a  more  continuous  result  does  not  mean  either  solving 
the  present  equations  more  accurately  or  approximating  the  governing  equa- 
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Figure  21:  Analytical  approximations  for  the  transmission  coefficients  for 
e  =  0.08  and  rj  =  1%  at  large  k  (which  corresponds  to  large  angles  of  propa¬ 
gation).  We  have  plotted  the  coefficients  related  to  the  phase  integral  trans¬ 
mission,  exp(— tt/3+/c),  in  dark  type.  The  transmission  coefficients  found 
from  the  embedded  Budden  equation  for  tunneling  through  the  “gap”  are 
indicated  in  light  type. 


tions  in  a  better  manner.  Rather,  it  means  incorporating  more  physics  into 
the  boundary  conditions  of  the  model.  For  example,  we  have  solved  a  prob¬ 
lem  in  which  the  coefficients  correspond  to  transmitted  wave  flux.  However, 
when  these  waves  reach  ionospheric  altitudes,  they  interact  with  a  conducting 
ionosphere.  Wave  reflection,  transmission,  mode  coupling  and  absorption  at 
the  magnetosphere-ionosphere  boundary  are  at  present  a  vigorous  subject  of 
active  research.  If,  for  example,  the  waves  were  completely  reflected  at  that 
boundary,  then  a  transition  at  a  critical  k  does  not  appear  because  waves  are 
reflected  at  all  k  and,  hence,  there  is  no  discontinuous  jump.  If,  on  the  other 
hand,  the  transmitted  waves  are  absorbed,  then  we  would  expect  an  abrupt 
transition  at  some  value  of  k  which  corresponds  to  reflection  at  the  altitude 
of  the  absorbing  layer.  Moreover,  density  gradients  play  a  substantial  role  in 
determining  the  exact  location  of  the  reflection  layer  and  thus  of  the  exact 
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Figure  22:  Coefficients  for  large  k  with  e  =  0.08.  Strong  absorption  is 
peaked  about  r/  =  1. 


behavior  of  the  coefficients  with  k. 

As  an  example  of  the  importance  of  boundary  conditions,  suppose  we 
were  to  assume  that  waves  incident  on  the  ionosphere  are  completely  re¬ 
flected.  Then  in  terms  of  the  phase  integral  analysis  discussed  in  part  1, 
we  find  trivially  that  for  an  incident  upgoing  ion-cyclotron  wave  from  below 
the  ion-gyrofrequency  the  amount  of  wave  flux  Tjc  =  exp(  — 7 r/?+/e)  contin¬ 
ues  upward.  Because  this  mode  is  completely  subdominant,  it  couples  to  no 
other  modes  and  hence,  the  remaining  wave  power  is  completely  absorbed. 
Thus,  for  example,  in  the  case  of  Figure  12  where  no  absorption  occurs,  we 
would  find  that  60-70%  of  the  flux  transmitted  on  the  coupled  branch  as 
Cjfs  would  be  reflected  at  the  ionosphere  and  then  absorbed  at  the  oxygen 
resonance  so  that  in  actuality  the  coefficients  which  we  have  obtained  serve 
as  a  lower  bound. 

In  sum,  our  results  provide  insight  with  regard  to  a  lower  bound  estimate 
of  the  absorption  near  the  resonance-cutoff-crossover  frequency  triplet.  The 
results  will  be  somewhat  modified  if  we  include  more  physics  into  our  analysis, 
but  it  is  certainly  clear  that  substantial  wave  power  is  available  to  heat 
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minority  species  ions  such  as  oxygen. 


5  Summary 

In  order  to  incorporate  the  presence  of  coupling  both  above  and  below  the 
heavy  ion  gyrofrequency,  check  our  analytical  results,  and  extract  the  values 
of  all  coefficients,  we  have  solved  the  system  of  coupled  differential  equations 
(2.1)  numerically.  In  order  to  integrate  the  equations  around  the  singularity, 
we  incorporated  physical  considerations  into  our  analysis.  This  prescrip¬ 
tion  involved  continuing  the  multivalued  solution  to  the  differential  equation 
“above”  the  pole.  By  imposing  boundary  conditions  on  the  numerical  solu¬ 
tion  consistent  with  downcoming  waves  in  a  particular  mode,  we  were  able  to 
determine  the  coefficients  associated  with  wave  propagation  near  the  heavy 
ion  gyrofrequency. 

We  found  that  our  numerical  solutions  are  in  reasonable  agreement  with 
the  analytical  solutions  that  we  have  obtained.  Incident  magnetosonic  waves 
couple  strongly  near  the  crossover  frequency  for  reasonable  values  of  the 
minority  species  density,  and  for  large  angle  of  propagation  the  waves  are 
completely  reflected.  Strong  absorption  (the  order  of  10-20%)  occurs  over  a 
substantial  portion  of  the  spectrum  peaking  near  perpendicular  propagation. 
Ion-cyclotron  waves,  on  the  other  hand,  undergo  Budden-like  tunneling  at 
the  cutoff  frequency  for  nearly  parallel  propagation  and  at  the  ion-ion  hy¬ 
brid  frequency  for  nearly  perpendicular  propagation.  The  absorption  is  very 
strong  (~  20%)  over  most  angles  of  propagation  for  reasonable  values  of 
the  minority  species  concentration.  In  either  case  the  absorption  levels  are 
consistent  with  the  power  required  to  heat  ion  conics. 
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6  Appendix  A:  Details  of  Numerical  Calcu¬ 
lation 


In  the  asymptotic  region  where  the  WKB  approximation  is  valid,  the  solu¬ 
tions  to  Eq.  (2.1)  are  close  to  the  set  of  basis  functions 


0i  =  n^2  exp  -  J  n+(z)d 
02  =  n+1/,2exp  y  n+(z)d; 
03  =  nZ1/2  exp  -  J  n_(z)d: 
04  =  nl1^2  exp  ^+-  J  n-(z)di 


(2.A1) 


which  are  asymptotically  equivalent  to  the  WKB  solutions.  The  functions 
n±  are  solutions  to  the  dispersion  relations  associated  with  Eq.  (2.1) 


(2.A2) 


and  have  nontrivial  spatial  dependence  according  to  the  function,  g.  For 
parallel  propagation  these  WKB  solutions  correspond  to  a  downgoing  LHCP 
wave,  upgoing  LHCP  wave,  downgoing  RHCP  wave,  and  an  upgoing  RHCP 
wave  respectively.  Indeed,  these  functions  are  equivalent  to  the  simpler  func¬ 
tions  0|,  0“,  01  and  01  respectively  of  part  1.  This  set  of  WKB  solutions 
provides  a  suitable  set  of  basis  functions  in  the  asymptotic  regime  far  away 
from  the  pole  in  the  complex  plane. 

The  Poynting  flux  (2.10)  in  terms  of  the  upper  and  lower  components,  ip 
and  <p,  of  (2.1)  is 

J  =  i{ipip*'  +  (pep*'  -  ip*ip'  -  ftf)  (2. A3) 

In  the  region  where  the  WKB  solutions  adequately  describe  the  actual  so¬ 
lutions,  we  find  that  the  coefficients  of  the  WKB  solutions  <p  and  ip  are 
algebraicly  related  by  Eq.  (2.1).  Then  to  first  order  in  the  WKB  parameter 
the  Poynting  flux  above  is  proportional  to 

J+  =  -htft\2  +  It Ut\2  -  I73/-I2  +  It4+/-I2  (2.A4) 
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and  the  Poynting  flux  below  is  analogously 


J~  =  -l7f/;l2  +  It2"/+ I2  -  l73“/r|2  +  |t4-/:I2  (2.as) 

where  the  factors  /  describe  the  relationship  between  the  WKB  representa¬ 
tions  of  ip  and  <p.  That  is 


/±  =  1  +  (3(n±  ~  r  +  KV2)) 


(2.A6) 


and  the  superscript  j  refers  to  the  location  at  which  this  expression  is  to  be 
evaluated,  namely  in  the  asymptotic  region  where  the  n*  have  been  deter¬ 
mined. 

The  power  flow  can  be  identified  from  the  various  terms  in  Eq.  (2.A4, 
2.A5)  so  that  for  the  case  of  magnetosonic  incidence  as  in  Eq.  (2.11)  the 
coefficients  are 


Tms  — 

Rms  = 

W2 

Cms  - 

wn/ft  i 

Cms  — 

(2.A7) 


For  a  downcoming  ion-cyclotron  wave  as  in  Eq.  (2.12)  we  find 


Tic  = 

wjun  i 

Ric  = 

Kl2 

III 

O 

Cfc  = 

(2.A8) 
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